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SYNOPSIS 

An uppei — bound solution methodology has been developed for 
the cold densifi cation rolling of sintered porous metal strips. 
Based on the recently derived plasticity theory for porous 
metals, the density and velocity fields during rolling have been 
estimated by minimizing the rolling power. The effects of 
rolling parameters such as friction between the rolls and the 
strip, roll diameter, roll speed, strip thickness and percentage 
reduction and also the material characteristic such as initial 
relative density of the porous metal, on the densif ication 
behaviour and rolling power are investigated. The other process 
variables such as roll pressuure and percentage slip are also 
calculated. The multi-pass cold densif ication rolling with and 
without intermediate annealing is analysed. The effect of work 
hardening is also taken into account. The density and strains 
have been verified experimentally for multi-pass rolling of poros 
coppef strip. 

In the upper bound analysis, a stream function form which 
produces a kinematically admissible velocity field is assumed. 
The total power for rolling, which is the sum of the internal 
power for plastic deformation, frictional power and shear power 
losses along the surface of velocity discontinuity, is minimized 
to find the velocity and density distributions. It is also shown 
that, there exists an optimal entrance strip velocity which 
corresponds to minimum power consumed. 



The distributions of density, velocity and strain indicate 
that the deformation is homogeneous. The results for roll 
pressure clearly show that the roll pressure is much smaller for 
porous metal than that of solid material under similar rolling 
conditions. It is seen that the exit relative density of strip 
increases with increase in friction factor. Also, the power 
required to deform the material is higher for large friction 
factors. The roll speed and roll radius do not affect the 
densi f i cat i on . However, the results show a marked increase in 
the power consumed at higher speeds. The effect of initial 
strip thickness is also investigated. It is shown that 
densifi cation is larger for smaller strip thickness. It is 
also observed that the densif ication behaviour is significantly 
influenced by initial relative density of strip. 

The calculated results for density and strains are in good 
agreement with the experimental ones for multi-pass rolling of 
porous copper strip. It is also shown that the densif ication of 
porous copper strip is not affected fay the mode of cold rolling 
i.e. cold rolling with intermediate annealing or cold rolling 
without annealing. However, the work almost increases by factor 
of two, if annealing is not performed. 
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CHAPTER I 
INTRODUCTION 

1.1 ROLLING OF POROUS METALS: 

For staking products via traditional powder metallurgy 
technique, powders are compacted in dies and subsequently 
sintered at elevated temperatures. The products thus obtained 
usually contain porosity of the order 3-10%, and are 
characterised by p'oor ductility value, which limit their usage to 
low stress applications. Various methods such as forging, 
extrusion and rolling are used to improve the strength, and in 
particular ductility, of the sintered product. Hot or cold 
rolling of sintered green strip is done to produce strip of full 
density. 

The entire process of making strip by powder metallurgy 
method consists of making a green strip from metal powder. This 
is generally done by roll compaction of loose metal powder. 
Alternatively green strip can also be obtained by slurry method 
C13. The green strip is then sintered at the appiropir iate 
temp'erature to increase the strength. Subsequently the sintered 
strip is densified by either hot rolling in one pass or repeated 
cold rol 1 ing— anneal ing cycles. Finally, the strip is cold rolled 
and annealed to improve surface finish and mechanical properties. 

The deformation behaviour of sintered porous metal is 
greatly influenced fay the presence of voids. The analytical 
methods that are used to study the deformation of ordinary pore— 
frefe metals cahndt be Applied directly to porous metals, due to 
the additioifal pV'o'b’lems'’ created by thfeae veids.' Kuhn C23 and 



Green C3I1 independently proposed yield criteria and stress-strain 
relations, which take into account density changes during the 
deformation of porous metals. However, the application of these 
theories was limited to uni-axial compression and tension or 
plane-strain compression, and could not be extended to analyse 
the complicated deformation processes. Shima and Oyane II4I] later 
proposed a plasticity theory by extending the yield criterion to 
a general stress state. 

Most of the studies carried out in the area of powder 
metallurgy processing of fully dense strip have been concerned 
with the feasibility of the process for producing a sound 
product. A number of analytical studies have also been carried 
out on roll compaction of loose powder to produce green strip. 
However, very little information is available on the 
densifi cation rolling of a sintered porous strip. The present 
study concerns with the cold densifi cation rolling of a porous 
copper strip and describes a theoretical methodology to analyse 
the various aspects of material deformation and densif ication 
characteristics . 

1.2 REVIEW OF THE PREVIOUS WORK: 

The densif i cation of sintered P/M products have been widely 
studied experimentally. However, only a few theoretical studies 
on this area are available in the literature. Most of the 
earlier analysis, have dealt with forging of sintered powder 
compacts. While, considerable amount of work has been done on 
solid metal strip rolling for a wide range of rolling parameters, 
very little information is available at present on the 



denei f ication rolling of poroua metal etripa. So far no attempt 
has been made to study the importance of various rolling 
parameters on the densifi cation rolling of P/M sintered strip. 

A major objective of mathematical modelling of metal forming 
processes is to provide the necessary information for proper 
design and control of these processes. Therefore, the method of 
analysis must be capable of determining the effects of various 
process parameters on the metal deformation characteristics. A 
wide spectrum of solution techniques have been emplouyed in the 
past to analyse the metal forming processes. The most popular 
among these are the upper bound analysis and the finite element 
method. The finite element method has received considerable 
attention as being the most powerful technique for the analysis 
of metal forming processes. The major advantage of the finite 
element method is that the method can be applied to a wide • class 
of boundary value problems without restrictions of workpiece 
geometry. It is possible to take into account the realistic 
boundary conditions and material properties. However, the amount 
of computational effort required for FEM analysis is often very 
high. The upper bound technique is also a useful analytical tool 
for metal forming problems. One of the major advantages of this 
method is that it is relatively simpler than FEM and yet is 
powerful enough to predict the values of important design 
parameters or operational variables. 

The rolling of fully dense material has been analysed 
extensively by various authors using the upper bound theorem or 
the finite element method C5-143. The effects of rolling 
parameters- such as friction tfactor^ roll diameter r roll speed. 
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percentage reduction, front and back tension etc. have been 
studied in detail. Avitzur C5J analysed the cold strip rolling 
of an ideal plastic material using the upper bound approach. The 
conditions of constant Columb friction coefficient were 
investigated. Recently, Avitzur and co-workers CfeD have proposed 
the use of two triangular velocity fields for the upper bound 
solution of strip rolling. The effects of various process 
parameters upon the total power has also been studied in this 
work . 

The finite element analysis of cold strip rolling was 
attempted by Rao and Kumar C7J. A plane strain problem was 
treated using the incremental displacement method. The stress 
and strain histories along the material path were traced and a 
nonlinear analysis of the problem was carried out. The 
deformation and pressure distributions in the plastic zone were 
obtained. The elasto-plasti c finite element modelling of cold 
strip rolling was investigated by Liu and co-workers C8,93. 
Assuming plane strain conditions, the elastic and plastic 
deformations both inside and outside the nominal contact zone 
have been predicted. Solutions for nonsteady and steady state 
rolling have been obtained for a wide range of rolling 
parameters. Mori et al. C10D have simulated the deformation 
characteristics for plane strain rolling by the rigid plastic 
finite element method. A slightly compressible material with and 
without work hardening has been considered. Both steady and 
nonsteady state strip rolling have been analysed. The three- 
dimensional deformation in rolling has been studied using FEM by 
Mori and Osakada C:113. A constant coefficient of friction have 
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been assumed between the plate and the rolls. 

As compared to solid metal rolling, the strip rolling of 
porous material has not yet been analysed in so much detail. 
Mori and Osakada C153 developed a finite element model for 
analysing the plastic deformation of sintered powder metals. The 
plane strain hot rolling of sintered powder plates were 
considered and the distribution of density, strain, strain rate 
and stress in the plate as well as roll pressure, were 
calculated. Using a variational formulation by rigid plastic 
material along with the plasticity theory of Shima and Oyane C4T, 
they performed a transient analysis in which the calculation 
was repeated from the initial state till the final steady state 
rolling was achieved. The non-steady state rolling was treated 
because the authors experienced difficulties in achieving a good 
convergence for the steady state scheme. It was found that, the 
computed stresses and strains were greatly influenced by the 
density distribution. In other related studies C16,17D the same 
authors have studied upsetting and the indenting of cylindrical 
billets. They have employed the rigid plastic finite element 
method to predict the stresses and the relative density 
distributions within the deformed material. 

It is clear from the review of literature presented above 
that considerable gap exists in the theoretical understanding of 
cold densif icatlon rolling of porous metal strips. In 
particular, simple analysis based on the upper bound technique 
are not available. The present study is a modest attempt to 
bridge the gap in the theoretical knowledge on this subject. 
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.3 OBJECTIVE AND SCOPE OF THE STUDY: 


< 3 ) 


' i i ) 


( i i i ) 


( iv) 


< V) 


(Vi ) 


The .present study aims to achieve the following objectives: 
To analyse the densifi cation behaviour of porous metals 
during cold rolling, using the upper bound technique. 

To provide estimates of the total power required for 
rolling and the break-down of rolling power into its 
various components. 

To predict the variations of velocity, strain 
components, and density within the deformation zone, 
and roll pressure and percentage slip. 

To conduct a detailed parametric study and investigate 
the effects of process and material parameters such as 
roll diameter, roll speed, friction factor, percentage 
reduction, strip thickness and the inlet strip density. 
To analyse single-pass and multi-pass cold rolling; 
studying multi-pass rolling with and without 
intermediate annealing. 

To verify the theoretical predictions on densif ication 
and strains experimentally. 


Due to limitations on the range of published data on 
material properties the scope of the study has been restricted as 
fol lows: 


(i> Simulation has been carried out for copper strips in 

the relative density range of 0.7 < P < 0.975. 

(ii> The range of parameter values considered are: 

roll radius: 67.5 mm - 350 mm 

roll speed! 550 mm/s — 800 mm/s 



friction factor: 


0.3 


0.6 


percentage reduction: single pass: 57. - 357. 

multi-pass : upto 86.67. 

( theoreti cal ) 

upto 61.4% 

( experimental ) 

initial strip thickness: 0.5 mm - 6.6 mm 

<iii) Steady state rolling has been analysed. The stress- 
strain relationship proposed by Shima and Oyane C4ZI have 
been considered. 

(iv) Density has been assumed to be uniform over each cross- 
section within the plastic zone. A power function form 
of density variation with the strip thickness has been 
assumed as proposed by Bhargava [1183. 

(v) Only densif i cat ion and strain characteristics have been 

' measured experimentally. 
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CHAPTER II 

MATHEHATICAL MODEL OF COLD DENSIFICATION ROLLING 

2.1 PLASTICITY THEORY FOR POROUS METALS: 

Powder metallurgical techniques, combined with conventional 
deformation processes, have become quite useful for the 
fabrication of various engineering products. Sintered powder 
metals are employed as starting materials in many of the common 
working processes such as forging and extrusion. In the 
modelling of these deformation processes, in general , plasti city 
theory has been found to be an invaluable tool for the 
calculation of forming loads and stress distributions. Powder 
metal working processes, however, do not come under the purview 
of common plasticity theories, since during the plastic 
deformation of porous metals, the material volume undergoes a 
drastic change. Moreover, conventional yield criteria are not 
applicable and the constitutive relations between stress aand 
strains are complex with a strain dependence on density. 
Analysis of powder metal working, therefore, requires the use of 
special plasticity theories. In recent years, the plastic 
deformation character ist i cs of porous materials have aroused 
wide research interest. Yield criteria and stress-strain 
relations for porous materials in uniaxial tension/compression 
have been proposed by Kuhn CSU and Green C33. These have, later 
been extended to a general stress state by Oyane C43. 



E.1.1 Yield behaviour of P/H materials: 

It is a well-known fact that plastic deformation of a fully 
dense material does not involve any volume change. Also, 
hydrostatic stress has no influence on the yield behaviour of a 
fully dense material. The yield criterion for such a material, 
is therefore, only a function of the second invariant (J’g) of 
the stress deviator throughout the expression 

f = (3 (2.1) 

where, f is the yield surface. The invariant J’g is, in turn 
given by the expression 

J ’2 ~ ~ (dg —dg)^ + (dg “ d-j)^Ii/6 

( 2 . 2 ) 

with d^ , Og and Og being the principal stresses. It is evident 
from the equation (2.2) that yielding depends only on the stress 
difference for a dense material and that a pure hydrostatic 
stress will not influence yielding. 

In the sintered powder metal, however, there is a change in 
material volume after deformation and the yield behaviour if 
sensitive to the hydrostatic stress. A criteron was proposed by 
Green C33, which was based on elastic-plastic deformation of a 
volume element containing spherical pores. The criterion is 
given by 

6 = J’g + a j2 (2.3) 

1 

where, a and 6 are the functions of the porosity and Y is the 
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yield strength of the material. However, it was found that the 
above relationship overestimates the Poisson’s ratio at low 
densities and predicts a much lower sensitivity to increasing 
density as compared to experimental measurements. The Poisson’s 
ratio for solid metal is equal to 0.5, while it is smaller for 
porous metals and increases till it attains a value of 0.5 when 
the density approaches 100X. 

Another criterion for the yielding of porous materials as 
proposed by Kuhn and Downey was based on the observed 
relationship between and density C20. This relationship is 
given by 

lateral strain 

y = = 0.5 P® (2.4) 

longitudinal strain 

where, a = 1.92 for room temperature deformation and a = 2.0 for 
hot deformation. Kuhn and Downey modified the yield criterion 
for a fully dense material and proposed the following criterion 
for porous materials! 

f = C3 J’s - (1 - 2y) Jg (2-5) 

2 

where, Jo is the second invariant of stress given by J’g - J /3. 

1 

The effect of hydrostatic stress was also included in the form 
of a quadratic term through Jg ®nd thus, satisfying the 
requirement that the yield criterion should be an even function 
of hydrostatic stress. Evidently the second term in equation 
(2.5) vanishes as v 0.5, for a fully dense material. 



Both the theories so far discussed have been used only for 
simple stress states such as uni-axial compression and tension or 
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plane-strain compression. It has not been possibale to apply 
these approaches for the analysis of practical deformation 
processes. 

Shima and Oyane C43 later proposed a yield criterion which 
has been successfully applied to analyse complicated deformation 
processes. The present theoretical simulation is also based on 
plasticity theory derived by Shima and Oyane. The yield 
criterion proposed by these authors is given by 

F Cf(ori + (0^ -ai)^>/2 + ( ff^/f 

( 2 . 6 ) 

where, the factor f represents the degree of influence of the 
hydrostatic stress component on the onset of yielding of the 
porous material and it may be a function of the relative density. 
F in the above equation can be related to the yield stress of the 
matrix metal5^q, which is a materiala parameter. The equation 
(2.6) can, therefore be rewritten as, 

faeq = CC(Cr^ - c ^)2 + (0^ -c^)^ + <(^3 ’ +( c^/f 

(2.7) 

In the above equation, f* represents the ratio of the apparent 
stress applied to the porous body and the effective stress 
applied to the matrix and is again a function of the relative 
density. The functions f and f’ can be determined experimentally 
C43 from uni-axial compression test. For porous copper. 
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(E.S) 


2.49 ( 1 - P 

and 

f = p2.5 (2.9) 

However, these forms of expressions for f and f* are also 
applicable to sintered iron and sintered alluminium. 

The yield surface of a porous body given by equation ( 2 . 7 ) 
is an ellipsoid, the major axis of which coincides with the c 
axis. On the other hand, for solid materials, the yield surface 
is a cyalinder whose axis coincides with axis. 


2.2.2 Stress-strain relations: 

For a porous material, the principal strain increments d e i . 
d Cg and d£^ can be obtained by partially differentiating the 
plastic potential (g) with respect to cr^ , and 03. Thus, 


39 5, 

dC^ = dX’ --- = dX C 0-, - (1 - 2/g f*^) 0^ > 

1 


3g P 

deg = dX* = dX C 0g - ( 1 - 2 /g f*^) 0^, > 


9 9 2 

deg = dX’ = dX foj - (1 - 2/g f^) 

The volumetric strain d y can be written as 


( 2 . 10 ) 


dp 2 

d£y = d + deg + dC3 = ~ ~o~ ~ dX( 2/3 f ) 0|n 


( 2 . 11 ) 



where, dX is a proportionality constant. If Wp is the 
incremental plastic work done per unit volume of the porous body, 
then 

dWp = + 02 dEg + de3 (2.1E) 
Alternately, 

dWp = ‘^^eq (E.13) 

where, d Ceq equivalent strain corresponding to the 
stress state of a^q. Substituting equation (2.10) into equation 
(2.13) and rearranging the terms, we have 


dX 


3 P 
2(f • 



(2. 14) 


Finally, d£gq can be written as 

f ’ 

de^q = t:(2/9) C (dei-d€2)2+(d€e-de3)^+(de3-dei)^> + (fdey)^3''^^ 

(2.15) 

Summarising the above equations we have 

- 1 

^eq = ""sTi" +(<^2 " <Ji )^>/2+(9;.,/f 


de-j = - 


2 P 


1 

2<2.5)-1 


dc, 


eq 


eq 


C 0f<l- 2/9f^) 


de 


de 


2 ' 


H P 


2(2.5)“1 


eq 


eq 


< ag - (i-2/9f ) a^y 


(2.16) 
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(E. 17) 


des = 


.S(E.5)-1 


®eq 




2/9f‘^) or^, > 


d £y = 


dH 

T~ 


de, 


eq 


,2(2.5)-1 ) 


m 

^2‘ 


eq 


and further 

dEgq = C(2/9)<:(de^-d€2)^ +(d€e-de3)^+(d Cg-d £ ^ )^> + ( f d 

<2.18) 

when the relative density ( P ) •* 1, then f , ft,1 and de^ = 0, 
and these equations reduce to the well-known stress-strain 
equations of conventional plasticity theory of a fully dense 
material . 

2.1.3 Stress-strain curve for the matrix material: 

In the basic equations of plasticity theory for porous 
materials given above, Ogq and 6gq refer to the equivalent 
stress and the equivalent strain respectively, to which the 
matrix material is subjected. Therefore, if a stress-strain 
curve for a sintered material is known, it can be converted to a 
stress-strain curve for the matrix material by utilizing the 
basic equations <2.16 - 2.18) C4I]. 



15 


The flow curve of the copper matrix was obtained from the 
simple compression test of sintered copper cylinders by Mori and 
Osakada II15D. The expression for the flow curve for copper is 

Cfeq = 431 ( E gq + 0.01)°'^^ MPa (2.19) 

The present theoretical simulation has been carried out in the 
relative density range 0.7 < P < 0.995. The plasticity theory 

proposed by Shima and Oyane is valid in this range of relative 
density. 

2.2 UPPER BOUND TECHNIQUE FOR POROUS METALS 

The upper bound technique is one of the simplest and most 
powerful tools used for analytical modelling of most of the metal 
forming problems. The upper bound theorem is based on a 
kinematically admissible velocity field which minimizes the total 
power required for the process. A kinematically admissible 
velocity field, in turn, must satisfy the following conditions: 

(i) The velocity components and their first derivative must 
be continuous except at allowable surfaces of velocity 
discontinuity. 

(ii) The velocity boundary conditions required by the geometry 
must be met. 

The upper bound solution provides a value for the power 
equal to or greater than the actual power required to deform the 
material. The conditions in the solution procedure are relaxed 
to the extent that the equations of equilibrium, the stress- 
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strain relations, the stress boundary conditions are not 
necessarily satisfied at every point. However, the total force 
applied by prescribed surface tractions is a satisfied parameter. 
There are infinite number of k inemat i cal ly admissible velocity 
fields that satisfy the velocity boundary conditions. The true 
solution is the one which corresponds to the minimum power. In 
the application of the upper bound solution, the plastic 
boundaries of some chosen forai are assumed which define the 
plastically deforming region. Once the plastic boundaries are 
fixed, a kinematically admissible velocity field is assumed which 
satisfies the velocity boundary constraints. In this regard, a 
physical understanding of the flow of metal is quite useful. The 
derivation of the upper bound theorem for the plastic 
deformation of a porous metal is briefly discussed below. 

The yield surface for a porous metal is an ellipsoid in the 
principal stress space, which is concave to the origin. Thus, 
the principle of maximum plastic work can be applied which states 
that the true increment in plastic work 6w required to bring 

1 

about a given set of strain increments dEjj is greater than the | 

c ^ ! 

increment in plastic work OW which would abe required to bring : 

I 

about the same strain increment by any other distribution of | 

I 

stresses throughout a volume V, which confirms to the same yield 
criterion. Thus, 

6w - 6w* = / (Ojj - i e.j jv 

In terms of strain rate and rate of work, [ 

W - W* = dEij dV 



The principle of waximuin plastic work leads to 


( 




* 

Oij ) d dV > 0 


(S.EO) 


This criterion is always satisfied by the rigid plastic 

material, if its yield locus is convex outwards. In equation 

(E.EO), is the exact stress or-the true stress condition 

* 

during the deformation, a^j is any other stress which does not 
violate the yield criterion and is the exact stra.in rate. 

Thusi we have 


/ 

V 


^ij 





dV 


(E.E1 ) 


The prescribed velocity boundary conditions and surface tractions 
for any plastically deforming region are shown in Fig. E.1. 
Using the equation (E.13), the above equation can be rewritten as 


/ 

V 


* 

a 

ij 


e. . 

ij 


dV < /cr.. E.. dV 

“V 


^ ^dgq Egq dV 

V 


(E.SE) 


On a porous body of volume V and total surface A, the 
stress, strain rate, and other relevant quantities are denoted as 
f ol lows : 


= prescribed velocity over the portion of surface Ay 
= surface traction prescribed over the surface area Ap 
^ij “ exact strain rate 
= exact stress state 

Vj = any kinematically admissible velocity field satisfying 
* 

Vj = on Ay 




Vj = Prescribed velocity over 
surface Ay 

Fi = Prescribed traction over 
surface Ap 
ffjj = Exact stress state 

€ij =Exact strain rate 


FIG. 2-1 DEFINATION OF TERMS USED FOR THE UPPER BOUND THEOREM 
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* 
1 j 

* 


= strain rate derived from Vj^ 


a 


13 


= any stress which does not violate the yield condition and 

. « 

is connected to e | j . 


The principle of virtual work leads to 

# •» 

/ ajj dV = / Fi Vi dA 

V A 


(2.23) 


The equation (2.23) can be rewritten as 


/ o. . £ . . 

V 


dV = / Fj Vi dA + / 

Ay Ap 


* 

V,- dA 


(2.24) 


Setting Vi = Vi on Ay, the above equation becomes. 


/ Cij &ij dV = / Fi Vi dA + / Fi Vi dA 


V 

Now, 


(2.25) 


* , 

a. . r . . 
ij iJ 


* # 

Pa e 

eq eq 


= P<7eq Eeq 


Also, 

. « 

J a e 
V iJ iJ 


dV 


/ a. , e . , dV = / Pa 


V 


ij ij 


V 


eq eq 


Using equation (2.25) in equation (2.26) 


. * 

. dV 


(2.26) 


/ Fi Vi dA + / Fi Vi dA < / Po; 

Ay Ap “ V 


eq ^ eq 


dV 


or , 

/ 

Ay 


Fi Vi dA < / POeq ^eq i ’^i 


V 


Af 


(2.27) 


Equation (2.27) defines the upper bound theorem for a porous 
metal. The right han4 side qf this equation gives an upper bound 
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for the work rate of the unknown surface traction action on Ay. 
If the frictional shear stress is specified at the roll 
workpiece interface Aj, the upper bound theorem for porous metals 
can be stated by the following equation: 


» * «• 

i Fi Vj dA < / P Cfg e g dV - Fi Vi dA u* dA 

Ay - V -I Ap 

(2.28) 

where, Zi u is the relative velocity at the roll-workpiece 
interface. The power W, derived from a kinematically velocity 
field is given by the right hand side of equation (2.28). Thus, 

, * * 

W = / P0ga SgQ dV + Au* dA - / Fi Vj dA (2.29) 

V A-j- Ap 

The equation above describes the generalized upper bound 

theorem for a porous metal. It gives an upper bound for the 

total power required to plastically deform a porous metal. The 

true density and velocity distributions are obtained by 

minimizing this total power. 

2.3 MODEL ASSUMPTIONS: 

The following assumptions have been made in the present work 
while obtaining upper bound solution for cold densif ication 
rolling of sintered metal strips. 

(i) It is assumed that, the material remains perfectly 

rigid before it coriies under the foils and also after 
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emercee out. In other words, the deformation cone lies 
within the entry and exit sections as shown in Fig. 
(2.3) . 

(ii) A velocity discontinuity has been assumed only at the 
entrance due to the sudden change in the direction of 
metal flow lines. At the exit from the rolls, metal 
flow velocity is continuous and lies in the direction 
of the tangent to the roll surface. 

(iii) At each cross-section the relative density is uniform 

in the thickness direction because of the small 
thickness value encountered in strip rolling. 

Practically speaking, there is a slight variation in 
density across a section, with higher density near the 
roll surface as compared to that in the interior. 
However, this small variation can be neglected without 
much loss in accuracy. 

(iv) The lateral spread in strip rolling is very small as 
the width of the strip is large compared to its 
thickness. Hence, strain in width direction is assumed 
to be zero and rolling becomes a case of plane strain 
rolling. 

(v) Isothermal conditions are assumed as appropriate for 
the cold rolling of thin strips. 

(vi) The rolls are assumed to be perfectly rigid. 

2.4 DERIVATION OF KINEMATICALLY ADMISSIBLE FLOW FIELD: 

The Figure 2.2 shows a unit volume of porous metal 

undergoing progressive deformation along the flow line. It can 



be seen that the relative density increases in the direction of 
rolling due to compression and contact growth. Referring to Fig. 
(E.3), 2 one-II represents the deformation zone and g is a 
normalized cross-stream coordinate which varies from the value of 
zero at the axis to unity at the roll surface. If the streamline 
connecting and i y2 is to have smooth transition at the 
exit, the slope of this streamline should be zero at the exit 
section. It has been assumed that there is a velocity 
discontinuity at the entrance to the roll which results in a 
sharp turn of the streamlines at that cross section. The effect 
of velocity discontinuity has been taken into account during the 
evaluation of work according to upper bound theorem. 

In the deformation zone, each 1= constant curve represents 
a streamline. In order to derive the streamline shapes, it is 
important to note the facts that the roll surface is circulasr 
with radius R and due to symmetrical rolls, the mid-plane of the 
strip is straight. Thus, continuous smooth transition of 
streamline shapes is possible if one assumes. Streamlines of 
circular form with centre lying along the exit section are 
assumed. The streamline i = 1 represents the strip surface, 
which is a circular with centre placed at the roll centre. As one 
proceeds towards the mid-plane of the strip, the centre of the 
circular streamline moves away from the strip surface. The 
streamline g = 0 is the mid-plane whose radius is infinity, see 
Fig. (H.3). 

The assumed streamline pattern satisfies all the velocity 
boundary conditions. For instance, at the roll surface the 
normal velocity V|^ = O, this comes from the definition of 
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C ollapsed poros ity 

\P^=^ 


111 IV 


FIG. 2-2 PROGRESSIVE DEFORMATION OF A UNIT VOLUME OF 
POROUS METAL SHOWING CONTACT GROWTH 



: Initial relative density ! 
p 2 density ■ 


FIG. 2-3 SCHEMATIC REPRESENTATION OF VELOCITY FIELD 
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streamline, as there is no metal flow across the streamline. At 
the mid-plane and exit sections vertical component of velocity v 
= 0, because streamlines are horizontal. Thus, the concept of a 
streamline not only facilitates satisfaction of velocity boundary 
conditions but also gives a clear idea about the shape of flow 
lines in the deformation zone. 

The equation of a circular streamline can be written as 

(x - L)^ + <y - yo)^ = r^ (2.30) 

where , 

L = projected length of the arc of contact in the flow 
direction 

r = radius of the streamline under consideration 

y^j = y— coordinate of the centre of streamline arc = r + iyg- 


It 

may be 

noted 

that the coordinates of the 

centre of 

the 

streamline are 

given 

by (L, y^,). 

From the Fig. 

(2.3), it 

is 

evident 

that 

at X = 

0^ y = &yi- 

Substituting 

this value 

in 

equation 

(2.30) 

yel ds 

the following 

expression for 

the radiuus 

r 


of the streamline, 

S (yi - yg) 

r = + (2.31) 

2 £ ( y .j — yg ) 2 

The mass flow rate per unit width of the strip can be 
expressed as follows! 

f = t P ^ u-j (2.32) 

where , 

= initial relative density of strip, 

U.J = enJE-rance velocity of , ]the strip. 



y-| = entry thickness of the strip. 

For a compressible flow, it is possibale to define a stream 
function such that 


pu = ( ) (2.33) 

ay 

where. 


u = longitudinal or horizontal velocity of the strip at any 
point in the domain, 

P = the alocal relative density of the strip 

y == stream function = mass flow rate between a given streamline 
and the mid— pilane of the strip 


aw 

PV = - ( ) (2.34) 

a K 

V = transverse or vertical velocity of the strip at any point 
in the domain. 

P^, u-j and y.^ are constants and can be considered as inlet 
parameters. However unlike P-j and y^, the inlet velocity u^ can 
not be chosen arbitrarily for any rolling situation, as one does 
not have any external control in selecting a particular value for 
this parameter. In the present work, u-j has been treated as an 
unknown parameter which is determined from the upper bound 
theorem. 

Thus, defining C = P-i u.| y-j, equation <2.3E) can be 

rewritten as, 

W = C ^ (2.3S) 

The above equation indicates that mass flow rate is directly 
proportional to the streamline coordinate ^ . Combining 



equations (2.33 - 2.35), we can deduce following equations: 


3? 

3 X 

c =- 

3 X 

P V . 

II 

1 

c = P 

ay 

u 

Differentiating 

the streamline 

aradius r partially with 

o 

OJ 

1 

-l2 

3^ ^ <yi - y2^ 

ay 

2 ( y -1 - y2 ^ 

3y 2 


(2.36) 


~ y2 — ^ * <y ~ yo^^ 

ay 3y 


ay 

(2.37) 

(2.38) 


3£ 

<y-yo>>^'' “ y2 — ^ 

ay 


<yi-y2> 3^ 

+ ■3(r+y-y^^) 

2^ (y^-y2) 2 3y 


(2.39) 


H 

3y 


<y - yn> 


-L (yi - yg) 

Cy 2 (y-yo> + (r + y - y^) + >3 

2£ <yi - y2> 2 


(2.40) 

Similarly, differentiating equastions (2.30) and (2.31) 
partially with respect to x, and sirripl ifying 


as 

3x 


(X - L) 


-L (yi-yo> 

c y2<y-yo^‘^^'" y ~ — i + 

2S^(yi-y2) 2 


(2.41) 



di 

Now, eubetitutinfi for 

3x 

(2.41) into equation (2.36), 
V are obtained: 


3^ 

and from equations (2.40) 

3y 

the following expressions for u 


and 

and 


y - ( r + | yg) 


-L (yi-ys) 

y 2 (y-( r+Syg) ) + (y- I yg) C = + 3- 

2^ <yi-y2) 2 

(2.42) 


^i^iyi 


y£(y-( r+^y2> > + 


-L^ 

(y-Syg) 

2 r (yi-yz) 


j 

<yi-yz^ 

j 

2 


(2.43) 


It may be noted that at the exit all the streamlines become 
horizontal and this condition is satisfied as seen from equation 

(2.43) which gives V = 0 at x = L. The centre line is a 
streamline with ^ = 0 and the velocity component normal to it 
becomes zero because r •* » for 1=0. It can also be verified 
that for 1=1 (roll surface) the normal avelocity is zero. Thus, 
it is clear that the above velocity pattern represents a 
kinematically admissible velocity field. It may be noted, 
however, that the proper density distribution is not known as yet 
and it is determined by work minimization as stated by the upper 


bound theorem 
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2.5 NATURE OF DENSITY VARIATION IN POROUS METAL STRIP ROLLING: 

II has been shown by Bhargava C18D that the densifi cation 
rolling of powder metals is dominated by two processes, namely, 
the compression of particles which increases the 
inter parti cular contact area and the sliding of the p<articles 
at their neck regioans without causing any substantial plastic 

deformation. Considering homogeneous compression and grains with 
an initial spherical shape, he derived that the relative density 
variation can be represented by the expression 

1 - P 

= ( (2.44) 

1-^1 
c-_ > 

1 

where, 

P = local relative density 
= inalet relative density 
y = local strip thickness 
and , 

y^ = inlet strip thickness. 

The exponent b takes a value of 2 when there is no slip 
between the particles for the forming processes such as forging 
and upsetting. In the present work also, density is taken as a 
constant over each cross-section (homogeneous compression) with a 
power function form of dependence on the local strip thickness. 
The exponent b, however, is left as an unknown to be determined 
as a part of the solution. 

The assumption of uniform density over each cross-section is 
reasonably valid for the rolling of thin strips as used in the 
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present study. The value of b, in general, can be dependent on 
rolling parameters friction factor, percentage reduction, roll 
speed, roll diameter and initial strip thickness and relative 
density. For each set of parameters, the true value of b will be 
determined from the minimum work principle according to the upper 
bound theorem. This, in turn, will result in the true density 
distribution for the given set of rolling conditions. 


S.6 STRAIN-RATE CALCULATIONS 

The strain rates in the deformation zone are given by the 
following equations! 


3 u 

= strain rate in x-direction 

3 X 

3 V 


= = strain rate in y-direction 

3y 

• 3u 3 V 

~ + = shear strain rate (S.45) 

3y 3x 

The other strain rates are zero due to plane-strain rolling 

• m * 

condition, i.e. "iyz " ®z =0- It is very cumbersome to 

analytically -evaluate expressions for the strain rate components 
from the equations for u and v. Therefore, the strain rates are 
evaluated numerically in the present study. 


2.7 EVALUATION OF THE TOTAL ENERGY DISSIPATION RATE: 

The total power required to perform the strip rolling 
process is made up of three separate contributions. They are, 

(i) power expended for internal plastic deformation within the 



strip 
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ii) power to overcome friction at the strip surface. 

iii) power losses due to shearing along the surface of velocity 
discontinuity. 


Internal Power of deformation: 


The equivalent strain rate for porouus metal plane strain 
oiling is given by the expression: 


sq 


= P' 




-ey)^+(ey-e2) +(e2-Ej^) +3 Yx'y> + <f e 


(2.46) 


■iere, it, is the volumetric strain rate, is expressed as 




(2.47) 


id C, = 0, for the present case of plane strain rolling of 

irous copF'er strip also, for copper, we have 

1 

2.49 (1 - P 

ms, equation (2.46) becomes 


(2.48) 


>q 



:e equivalent stress can be written from the pilasticity theory 



•I u> 


'eq 



The power necessary for plastic deformation Up is given by 


3l 


Wp = IIP Cgq dx dy 


(2.51) 


The present ' analysis is carried out for porous copper in the 
relative density range 0.7 < P < 0.995. The equivalent stress- 
strain curve for copper C153 is expressed as 


**eq = *^31 < Eeq 0.01)°*^^ 

Using equation (2.52) in equation (2.52), 


MPa 


(2.52) 


Wp = 431 IIP 


“ a 9? 

Eeq < Ceq + 0.01)^’'^^ dx dy 


•eq 


(2.53) 

undergone by a 


The cumulative equivalent strain 
particle is determined by integrating the equivalent strain rate 
gq along the. streamliane corresponding to the motion of the 
particle. 


Friction Powers 

In porous metal rolling, the actual friction conditions are 
quite complex. The frictional condition is usually represented 
by a coefficient of friction or a friction factor. The 
frictional shear stress is a function of the contact pressure. 



location. Blip velocity etc. It reverses its direction at the 
neutral point. The frictional shear stress is usually expressed 
by the product of contact pressure and the coefficient of 
friction. The other simple approach for representing the 
interfacial friction phenomenon is through the prescription of a 
friction factor. The friction factor is defined as the ratio of 
the frictional to the yield shear strength of the material. 
Metal rolling is a friction aided process and therefore, one can 
expect a high value of friction factor. The value of friction 
factor can be chosen on the basis of the surface roughness of 
rolls and the lubricating conditions. The higher the friction, 
the larger the power needed to perform the process. The larger 
value of friction factor makes the rolling process more stable 
CS'T. However, it becomes less efficient from the point of power 
required for rolling. 

In the present work, the frictional resistance on the 
interface is assumed to be constant for the sake of simplicity, a 
constant friction factor is considered over the whole arc of 
contact. Also it has been assumed that both the rolls have same 
surface conditions. 

The frictional stress is expressed by 

T=mK <S.54) 

where, K is the criatical shear stress of the material whose 
value is obtained by setting the normal stress components in 
equation (2.50) to zero. Thus, 
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K = 


a p2-5 

eq ^ 


/■3 


(E.55) 


The velocity of slip at the roll-strip interface is. 


Ay = (U 




£ s . 
V ) ! 


i =1 


(2.56) 


where U is the roll speed. The frictional power is calculated as 


Wf = /tAV dl 


(2.57) 


The friction power is evaluated along the arc of contact. 
Substituting the relative velocity V into equation (2.57), we 
get , 

W. = / p2.5 cf <u -/u^ + v^) dl (2.58) 

/ 3 ^ 

assuming constant m along the arc of contact. 


Shear power losses at entry: 

A velocity discontinuity exists at the entry because all the 
streamlines bend suddenly at entrance to the roll. This 
discontinuity has been considered taking into account the 
straight boundaries of deformation zone. It has been assumed 
here that, the plastic deformation starts immediately after the 
strip enters the roll. In actual practice, plastic deformation 
starts even before the strip is gripped between the rolls and 
plastic boundaries are not straight. 

The shear losses at the discontinuity are given by 


W 


s 


1 

/ 

>^3 


“'eq 


2.5 


Av dl 


(2.59) 



where Lv is the difference in v component of the velocity at the 
discontinuity. 

Expression for total power: 

The total power required for rolling is given by, 

W = Wp + *W^ + (S.60> 

For a given reduction and friction conditions, the work is 
optimized with respect to exponent b and the entrance velocity u.-| 
of the strip. 



CHAPTER III 


NUMERICAL SOLUTION PROCEDURE 

3.1 OVERALL SOLUTION STRATEGY 

In order lo obtain the true density and velocity 
distributions through the upper bound technique, the primary 
quantity which needs to be evaluated is the total rolling power. 
This has been accomplished by numerical integration of the 
various components of power over the def oriTiation zone, using a 
kinematically admissible velocity field and a guess density 
distributions, given by equations (2.42 - 2.44). For a given set 
of rolling conditions, the minimum rolling power has been 
evaluated by systemati cal ly varying the two free parameters of 
the problem, namely: the exponent (b) of the density 
distribution and the strip velocity (u^) at the inlet cross- 
section of the deformation zone. The velocity and density 
distributions corresponding to minimum power have been taken as 
the true distributions. 

Due to the strong non-linear nature of the problem, it is 
not possible to evaluate the total power integral analytically in 
an exact manner. For this reason, numerical integration based on 
Gauss-Legendre quadrature has been resorted to. The numerical 
quadrature requires that for evaluating the area integral of 
internal plastic deformation power, the deformation zone be 
discretized into several small area elements? similarly, for 
calculating the friction and shear loss power contributions on 
boundary curves, these boundaries need to be discretized into 
several line segments. Both these requirements are met by 



discretizing the whole plastic deformation region into a finite 
element type mesh consisting of several small elements. All the 
elements have quadrilateral shape with 4 nodes each (see Fig. 
3.1). The nodes (which are placed at the vertices of the 
quadrilateral) aid in the computation of both the area and line 
integrals via numerical quadrature. A mesh consisting of S1 
elements in X-direction and 10 elements in y-direction (with 21 >: 
10 total elements) has been used in the present study for 
performing all the computations. Only half the strip thickness 
(upto the mid-plane has been considered due to symmetry). The 
shapes of the elements are such that two opposite sides of each 
element are aligned with particle stream lines while the other 
two sides are normal to the metal flow direction (See, Fig. 3.2). 
The total rolling power has been computed as the sum of the 
following quantities: 

(i) Internal plastic power over the 21x10 area elements. 

(ii) Friction power over the 21 line segments lying on the 
strip surface. 

(iii) Shear losses over the 10 line segments of the entry 
section . 

3.E CALCULATION OF FIELD VARIABLES AT NODES 

In order to evaluate the elemental integrals, it is 
necessary to find the values of the field variables density, x 
and y-direction velocities, strain components etc.) at the nodes 
of each element. This is easily achieved, if one recognizes the 
fact that ali the nodes lie on one streamline or the other. 
Since each stream line; is identified by a constant value of | 

I , i ,, , A . ■, 



4 


t>/ 


y 



X ,y =6lobal coordinates 
'n>^ = Local coordinates 


X 


G.3-1 ORDER OF GAUSS POINT SAMPLING POSITION FOR 
FOUR-NODED ISOPARAMETRIC ELEMENT 
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( stream-l ine coordinate), the field variables are sequentially 
computed as described below: 

(i) For the E1x10 mesh, the x-direction and ^-direction 

increments are respectively given by x = L/21 and 

AS = 1/10. From this, the x and S values at each node 

are calculated. 

(ii) Knowing x and S , from equations (2.5C>) and (2.3^!), the 
values of y and r for each node are determined. 

(iii) Using equations (2.42 - 2.44), the velocity components u 


( iv) 


and V and the density p are evaluated at all the nodes. 

* « 

In order to compute the strain rate components e^,, 
and Y j^y, and the cor responding cumulative strains 
€y and Y ;^y under gone by a particle, the following 


procedure 

is adopted. 

The strain 

rate 

3u 

components 
3v 3u 

related to 
3v 

the velocity 

derivatives 

7 

3 X 

3y 3 y 


— through the expressions of eq.(2.45>« 
3 X 


are 

and 


The velocity derivatives, in turn, are evaluated using 4- 
noded inter polation of u and v variation within each 

element. After obtaining the strain rate components 

• • 

Ey and Y ^y, the nodal values of cumulative strain 
components are evaluated by tracing the strain history of 
each particle along its stream line path. The 
incremental value of a particular strain component is 
obtained by integrating the corresponding strain rate 
component with respect to time. The evaluation procedure 
for the strain rate and strain components at each node is 
described in the next section in greater detail. The 



various steps involved in the evaluation of the true 
distributions are graphically shown in Fig. (3.3). 

3.3 CALCULATION OF STRAIN RATE AND STRAIN COMPONENTS 

FroiTi the nodal values of u,v and P , it is possible to 
obtain the values of these variables at any interior point of an 
elerrient through interpolation. For 4— noded quadrilateral 
elements, one can employ the following interpolation schemes 
within each element: 

4 

u = 2 Ni Ui 

i = 1 

4 

V = r Ni v^ (3.1) 

i = 1 

4 

P = 2 Ni Pi 

i = 1 

where Ui, Vi and Pj are the nodal values of the velocity 
components and density respectively, and Ni are the interpolation 
functions. It is more convenient to represent the interpolation 
fnctions in terms of normalized local coordinates ? and T] , 
rather than in terms of actual coordinates x and y. The 
expressions for the interpolation functions in local coordinates 
( iT , t}) are : 


1 1 

= (1-S) (1-Tl)s Ng (1 +?) (1 -tl) 

4 4 
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1 1 

Mg = (1+ ? ) ( 1+n ) ? N 4 ( 1 - ? ) ( 1 + T) ) 

4 4 

(3.2) 

The local coordinates are defined in such a way that they vary 
from -1 to +1 for all elements; thus, all quadrilateral elements 
(regardless of their shape) can be represented as squares of size 
2x2, in terms of their respective local coordinates. Such a 
feature is extremely useful in performing numerical integration. 

Now, for a quadrilateral element, it is possible to 

interpolate the global coordinates x and y at any point in terms 
of the values of nodal coordinates. 

It can be shown that 

4 

X = 2 l^i >«i 

i=1 (3.3) 

4 

y = S N^ y^ 

i = 1 

For the four-noded quadrilateral element, the conversion between 
the (x,y) and ( ? , T) ) coordinates can be achieved as follows. 
From chain rule , 
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The derivatives of interpolation functions aNj/a^ snd aNj/at) , 
in turn, are obtained by differentiating the expressions for Nj: 
in eq. <3.E) with respect to ? or T) . 

Now, the strain rate components can be evaluated using 
equations (3.1) - (3.8), as follows* 
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Lth the help of equations (3.4) - (3,11), it is now possible to 
Iculate the strain rates at each of the node points of an 
lement. However, an examination of the mesh geometry (Fig. 3.2) 
ows that most of the nodes are common to two or more elements, 
king into account this aspect, the strain rate components at 
ch node within the deformation zone have been calculated as 
.e appropriate average values between neighbouring elements. 



After evaluating the strain rate components the 
.orres ponding cumulative strain components are obtained as 
ollows. The incremental strains occuurring between two 
leighbouring nodes n and n+l lying on a particular stream line, 
ires 
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• • • 

here, ^x,av' ^y,av "Txy,av average values of the 
train rate components between nodes n and n+1, At is the time 
nterval required for a particle to travel from n to n+1, u^y is 
he average x-direction velocity between the two nodes and Ax 
s the step-size of the numerical mesh in x-direction. Using the 
bove expressions for incremental strain components, it is very 
imple to calculate the cumulative strains (or in other words, 
te strain-history) of material particles. Setting, e = 0, 
' y == 0 and Y j^y = 


0 at the inlet boundary (neglecting elastic 



or plastic strains the particles might have undergone before 
reaching this surface), the cumulative strains are progressively 
calculated along each streamline of the numerical mesh. The 
values of the strain components corresponding to all the nodes 
are stored as suitable arrays, for evaluating the stresses and 
the power integrals. 


3.4 ROLLING POWER CALCULATIONS 

As discussed earlier, the total rolling power is made up of 
three contributions, namely, the power for internal plastic 
deformation, power to overcome friction at the strip surface and 
the shear losses at the inlet boundary. The first one is an area 
integral over the deformation zone and the latter two are line 
integral computed on appropriate boundary curves. 

In the present study, Gaussian quadrature is employed for 
calculating both the area and line integrals. The application of 
Gaussian quadrature for computing area and line integrals is 
described below: 
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where, i is an index for the sampling points (also known as Gauss 
points) Fi orF(Xi, yj ) are the values of the integrand at the 
Gauss points and Ij or Ai represent the weighted length or area 
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corresponding to the i sampling point. In terms of the local 
coordinates ? and Tj , the evaluation of the integrals becomes 


even 

simple r . 

Thus , 





// 

(e) 

F dx dy 

+ 1 

= / 

-1 

+ 1 
/ 

-1 

F !J! 

m m 

dr dTj=: Z S WjWj 

i=1 j=1 

fF.!J!}. . 

* r J 







(3. 16) 

and 







J I 

( e ) 

F dl = 

L® 

+ 1 
/ 

-1 

F dr 
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= L® 2 Fi Wi (3.17) 

i==1 


where, Wj and wj are standard weights and L® is the length of a 
line element. The subscript is used in the above integrals 
indicates that the integration process is carried out over a 
quadrilateral element or the boundary line segment of an element. 
In the present study, consistent with the interpolation 
schemes employed, 2 sampling points per direction have been used. 
Thus, for the area integrals (2x2) areas points are used while 
line integrals are computed with 2 Gauss points, for each area 
element or line segment respectively. The weights and sampling 
point positions are summarized in Table 3.1. 
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Table 3.1 


i 

Wi 

i' n 

1 ' 1 

1 

1.0 

-1//3 

2 

1.0 

+ 1/^3 


Each of the power contributions can now be evaluated as 
follows. For the plastic power calculations, the dens i ty , 
equivalent stress, equivalent strain and equivalent strain rates 
are required at the Gauss points of each element. These are 
calculated using equations (2.45) - (2.52) and (3.9) - (3.14), 
along with interpolations of the form 


^ ^'^i X , i ’ ^ *^i i ' ^ xy 2 Nj Y xy , i 

i=1 i=1 i=1 


^ = 2 Nj ! Ey = Z Ni Cy,; ! Y 

i=1 1=1 


xy 


= 2 Nj Txy,l 
i=1 


(3. IS) 


applied at the Gauss points of the element. 


In the above expressions, i denotes the index for the nodes 
of an element. In order to facilitate the transformation of the 
area integral from (x,y) coordinates to ( ? ,t)) coordinates, it 
is necessary to calculate the determinant of the Jacobian matrix 
!J! at the Gauss points also. . 

Summing up the values of the product P*ffeq * ^eq the 
four Gauss points and multiplying with appropriate Gauss point 
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weights as shown in eq. (3.16), leads to the plastic deformatioan 
power contribution of one area element. Adding the contributions 
from all the 210 area elements provides the total plastic power 
over the whole deformation zone. 

The calculation of friction power at the roll surface 
involves the following steps. First, the frictional stress 
values and the local velocityship are determined at the Gauss 
points (which lie on the boundary) using equations (2.56) - 


(2.58) . 

The 

product L T is obtained at each 

of 

the 

two 

Gauss 

points 

and the frictional power contribution 

of 

one 

line 

segment 

is then evaluated via eq. (3.17). The sum 

of 

all 

such 


contributions from the 21 line segments on the roll surface, 
provides the total friction power. 

The shear losses at entry section are computed from eq. 
(2.59). The integration is performed in a manner similar to that 
of the friction power integral. The shear loss contribution for 
each of the 10 line segments along the entry cross-section is 
calculated and summed up, to obtain the total shear loss. The 
rolling power is then obtained as the sum of all the three 
individual powers. 

3.5 MINIMIZATION OF POWER 

The integral for the total rolling power W,has been 
evaluated numerically by Gass-Legendre quadrature, assuming an 
initial guuess for the value of b and u^. The derivative dW/db 
is evaluated by slightly perturbing the value of b by ^fa and 
using the central difference expression 



• • • 

dW W(b + ^b) - W<b - ^b) 

db Eiib 

To find the optitrium value of b, wesets 

dW 

= f (b) 
db 

and calculate the root of the expression: 


(3.19) 


(3. SO) 


f(b) = 0 (3. SI) 

by the Newton-Raphson technique. An improved guess value of b is 
obtained by the following equation: 
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where , 


f(b„) = 


(3. S3) 


df d^W 

db db^ 

In the above equatioOr the subscript n denotes the value 

t h * 

corresponding to n iteration. The second derivative of W in 

eq. <3. S3) is calculated through the central difference 

expression : 


S • 

d W 


db‘ 


W (b+^b) + W(b~ ^b) ~ SW (b) 


&b‘ 


(3.24) 


During the numerical implementation, . eq. (3.24) was converged to 
a value smaller than that of a prescribed parameter. Thus, 
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! f ( b ) ! <6 (3.25) 

was used as the convergence criterion for the solution, with 
chosen as equal to 0.0001. 

The entrance strip velocity u^, is then varied 
systematically to identify the, global minimum power for each set 
of rolling conditions. The distribution of density and velocity 
corresponding to the global minimum is then obtained and the 
corresponding results are taken as the true results for the 
given set of rolling conditions. The numerical computations have 
been performed over a wide range of rolling conditions, for both 
single pass and multi-pass rolling. For multi-pass rolling, the 
assuumption of no initial strains has been relaxed when 
intermediate annealing is not considered. For the calculations 
cor resp'onding to multi-pass rolling with intermediate annealing, 
the initial strains have been taken as zero and it is also 
assumed that annealing does not cause any densif i cation . The 
results predicted from the analysis are presented and discussed 
in Chapter V. 



r 



















53 


CHAPTER IV 

EXPERIMENTAL PROCEDURE 

Cold rolling experiments have been performed on stripis 
prepared from atomised copper powder. The powder sp>ecif ications 
ar e J 

Manufacturer - Greenback Industries, U.S.A. 

Copper content - 99.51% 

Powder size - 63.9% of -3E5 mesh 

Apparent density - 2.83 kg/m 

Hall flow meter rating - 3.45/50 gms 

The experiments have been carried out with copper strip having 
relative density (^) in the range 0.76 < P < 0.95. The copper 
strip of required density was initially produced by the die 
compaction technique. Approximately 140 gms of copper powder was 
taken and was mixed with 1.4 gms of methyl cellulose binder. 
The mixture was then compacted in a die of size 4.7 cm x 7.9 cm. 

The die and punch were lubricated with zinc stearate. The 

2 

compaction pressure used was 1451.65 kg/cm (load = 53.9 tons). 
The strip of dimensions 4.55 cm x 6.7 cm x 0.675 cm with a 
relative density of 0.7465 was obtained. 

The green strip thus produced was then sintered at a 
temperature of approximately 9E5°C in hydrogen atmosphere for a 
period of 70 minutes and then cooled in hydrogen for about 20 
minutes. After cooling, the weight and the dimensions of the 
strip were measured. It was found that there is some loss in the 
weight due to the removal of oxide film on the powder surface. 
The relative density of the strip after sintering was 0.7614 and 
the strip dimensions were 4.55 cm x 6.7 cm x 0.66E cm. 
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The sintered copper strip was cold rolled in a rolling mill 
with the following specifications: 

Roll radius = 67.5 mm 

Roll speed = SOS. 68 mm/s 

Horse power = 7.8 HP 

The rolls were not lubricated. The required amount of 
thickness reduction was carried out by multi-pass rolling. After 
each pass, density and dimensions of the strips were measured. 
Cracks were observed, after 30% reduction with respect to 
original thickness. The total reduction of 61.4% was carried out 
in seven passes. The true length and thickness strains were 
calculated from the obtained dimensions. The results for multi- 
pass rolling were summarized in Table 5*15. 
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CHAPTER V 

RESULTS AND DISCUSSIONS 

The results obtained through the upper bound technique for 
cold densif ication rolling of sintered porous copper strips are 
presented in this section. The main aspects studied here include 
the densif i cation behaviour of the porous strip and the power 
consumed during the rolling process. However, some effort has 
also been made to gain a fundamental understanding of the 
deformation process by analysing the velocity field, strain 
field, percentage slip and the distributions of roll pressure and 
relative density. The effects of parameters such as friction 
factor, percentage reduction, roll speed, roll radius and initial 
density and initial strip thickness have been highlighted. 
Results for multipass rolling with and without intermediate 
annealing have also been presented. 

Some of the calculated results have been verified 
experimentally. All the theoretical calculations have been 
performed for a sintered copper strip with relative density (P ) 
in the range 0.7 < P < 0.995, using the correlation for 
equivalent stress proposed by Shima and Oyane C4I]. 


5.1 IDENTIFICATION OF MINIMUM POWER CONDITION: 

As stated by the upper bound theorem (Chapter III), of all 
the kinematically admissible velocity fields, the true velocity 
field minimizes the total power. In the present study, in 
addition to the kinematically admissible velocity field, a 
density distribution witb an unkn.qwn axponent (b) is assumed. 



Thife exponent is determined by applying the minimum power 
criterion. 

Another important variable whose value is unknown a priori, 
is the strip velocity at the entry (u-j). A significant finding 
of the present study is that there exists an optimum entrance 
strip velocity which leads to the global minimum of power 
consumed for a given set of rolling conditions (Fig. 5.1). In 
view of the upper bound theorem, such an optimum entry velocity 
may be taken to correspond to the true velocity and density 
distributions. Therefore, for each set of rolling conditions, 
the exponent (fa) and the entrance strip velocity (u^) have been 
systematically varied through a Neweton-Raphson search procedure 
for identifying the global minimum power. The typical procedure 
employed to track-down the point of global minimum is illustrated 
in Fig. and Table 3.t. 

A close examination of the power contributions in Table 5.1 
shows that for u^ > u^^ and b = ^opt' plastic deformation 
power is high due to high strain rates; on the other hand, for u^ 
< u^ with b = bjjp^ friction power increases due to a high 
value of backward slip. When these two effects exactly balance 
each other, the optimal flow rate occurs (i.e. u-j = u-^^op^). 

It was generally observed that as the entrance strip 
velocity is varied, a global minimum in rolling power does not 
occur for reductions greater than 35%. It appears that for a 
fixed friction factor, the densif ication is improper at high 
percentage reduction (S5% and above) and the slip increases; 
these phenomena are probably responsible for no minimum existing 

However, the real reasons for the non- 


at high reductions. 



27.0 


57 



FIG, 5-1 THE GLOBAL MINIMUM OF POWER CONSUMED FOR A GIVEN SET OF 
ROLLING CONDITIONS 



existence of a global minimurri in our predictions for these 
conditions, are not clearly understood at this juncture. 
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5.E DISTRIBUTIONS OF THE FIELD VARIABLES IN DEFORMATION ZONE: 

The streamline patterns for two reductions (H07. and 57.) are 
shown in Figure 5.2. It is observed that streamline which are 
equally spaced at entry retnain equally spaced throughout the 
deformation zone. Ths feature clearly illustrates that the 
rolling process is homogeneous on account of the stnall strip 
thickness. The figures also indicate that for larger reduction, 
the arc of contact is longer and the streamlines are more curved, 
especially near the entrance section. This in turn, leads to 
larger strains. 

The u and v velocity distributions are presented in Figures 
5.3 and 5.4 respectively. The u-velocity distribution has been 
plotted only in terrris of the increment over the entry velocity 
(u^), for the sake of clarity. It is observed that the u- 
velocity is uniform while, the v-velocity decreases linearly from 
the surface to the centre, over any cross-section. The 
uniformity of u-velocity is a consequence of homogeneous 
compression. The v-velocity decrease towards the centre is a 
result of less bending of streamlines in the mid portion of the 
strip. In the rolling direction (x), u-velocity increases 
monot on i cal ly due to longitudinal elongation of material 
particles. The rate of elongation is high near the entry section 
where area reduction is also high. Towards the exit, the change 
in u-velocity becomes neligible. The v-velocity, on the other 
hand, decreases in the rolling direction. This is obviously a 
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Initial strip thickness =2 0mm 
Initial relative density i 0-70 
Roll radius =67-5nnm 




FIG. 5-2 DEFORMATION ZONE SHOWING STREAMLINES 
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FIG.5-3 U-VELOCITY DISTRIBUTIONS WITHIN THE DEFORMATION 
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Reduction s 20 V. Initial strip thickness =2-0mm 

scale s-10mm /s Initial relative density = 0-70 

Roll radius =67-5mm 

Roll speed = 508-68 mm /s 
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V-VELOCITY DlSTRBUTiONS WITHIN THE DEFORMA- 
TION ZONE 



FIG. 5 -4 


i^&ult of the streamlines becoming straight near the exit. At 
ower reduction, the distributions of u and v are very similasr 
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0 those at high reductions; however, the magnitudes of the 

hanges in u and v velocities between the inlet and exit are 

mailer. 

The strain distributions ( Ey and plotted in 

figures 5.5, 5.6 and 5.7. The normal strains in both 

longitudinal and lateral directions are uniform over each cross 
section. This feature also supports the conclusion of 
homogeneous compression, which was noted while discussing the 
streamline patterns. Both the normal strains increase rapidly 
near the inlet section and become almost constant towards the 
exit. This trend is clearly a consequence of the rapid changes 
in the velocity field near the inlet section. The shear strain 
is maximum at the surface and zero at the mid-plane, which is 
easily explained in terms of the larger bending of the 
streamlines near the surface. As expected, all the cumulative 

strains ( ' ^y undergone by a material particle, 

increase monotoni cal ly from the inlet to exit. The same strain 
behaviour is seen at all reductions, except that the magnitude of 

strains is larger at higher reductions. 

The variation of density, roll pressure and slip in the 
foiling direction are shown in Figures 5.8a and 5.Sb for the 
reductions of EOX and SX. The locations of the neutral point and 
the mauimuis roll pressure are also indicated in these figures. 
It is observed that the maKiEousi roU pressure occurs. close to the 
neutral point as expected. Although theory of rolling predicts 
that the maximum pressure should occur at the neutral point (due 
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Percentage reductions 20V. Initial strip thickness =2 0 mm 

Scale: 0-122 Initial relative density = 0-70 

Roll radius = 67-5 mm 
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FIG.5-5 DISTRIBUTIONS OF LONGITUDINAL STRAIN(Ex) WITHIN THE 
DEFORMATION ZONE 
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FIG.5-8a VARIATION OF RELATIVE DENSITY ROLL PRESSURE AND j 
PERCENTAGE SUP IN THE ROLLING DIRECTION FOR 20Vo j 
REDUCTION f 
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FIG. 5-8b VARIATION OF RELATIVE DENSITY, ROLL PRESSURE AND 
PERCENTAGE SLIP IN THE ROLLING DIRECTION FOR 
5 7o REDUCTION 
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to zero slip and zero shear stress), the upper bound analysis 
predicts a slight difference between the two locations. This may 
be due to the fact that stress balance is not applied everywhere 
within the deformation zone in the upper technique. It is seen 
that tl'ie variation of density is rapid near the entry and the 
curve becomes flat towards the exit. This trend is a direct 
consequence of the rate at which roll gap decreases with distance 
X, as reflected by equation (2.44) which governs the density 
variation. The variation of slip over the roll surface is such 
that the backward slip is larger in magnitude than the forward 
slip. The maximum backward slip which occurs at the entrance, is 
governed by the global minimum power criterion. It was also 
observed that the neutral point occurs closer to the exit than to 
the entry for the minimum power condition. With these facts in 
mind, since the velocity variation near the exit is slower, the : 

observation that forward slip is smaller than backward slip is 
easily understood. A comparison of the roll pressure curve for 
sintered porous copper with that of solid copper C93 indicates , 

that roll pressure is much larger for the solid material for a I 

given reduction and friction factor; this derives from the fact 

that the yield stress and consequently the work required for ^ 
deforming the solid metal, are higher than those for the porous 

material. At low reductions, the neutral point and the maximum 

pressure shift away from the exit section towards the mid-portion ^ 
of the strip. The magnitudes of both the backward slip and 

forward slip are smaller, this occurs chiefly because of the 

* * ♦ ETn*. ss niveh ihitial density/, smaller 

smaller arc of contact. For a given inn.i« 

4 ‘ '' ini t ^en'^ity. The ma>4imum 

percentage r^eduction leads to lower e«t density 
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roll pressure is also smaller which is contributed by factors 
such as smaller arc of contact, less dens i f i cat i on , lower strain 
hardening and smaller rolling power requirement in general. 

5.3 INFLUENCE OF ROLLING PARAMETERS 

(i) Friction Factor: The effects of friction upon the exit 
relative densdity are shown in Figure 5.9a. For the graphs shown 
in the figure, the initial relative density before rolling is 
taken to be equal to 0.7. For a given percentage reduction, the 
exit relative density of the strip is observed to increase with 
friction factor. This can be attributed to larger metal flow 
rate occurring because of better gripping conditions. When the 
friction factor is increased, the backward slip is reduced 
resulting in a higher metal flow rate. While, the forward slip 
exhibits a negligible variation (Figure 5.9b). However, the 
power required for rolling process increases with friction factor 
(Figure 5.9c). The increase in the rolling power is chiefly 
contributed by the frictional work. The results are shown in 
Tables 5.2 and 5.3. 

(ii) Percentage Reduction: Figures 5.9a,b,c and d also show 
the influence of percentage reduction upon densif ication, slip 
and rolling power. For any given friction factor, the exit 
relative density is larger at higher reduction. This is easily 
explained in terms of greater compression as well as the 
elongation of the particles in rolling direction. The elongation 
of the particles is an important mechanism of densif i cation for 
the high value of initial strip density assumed here 


It is 
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roll pressure is also smaller which is contributed by factors 
such as smaller arc of contact, less dens i f i cat i on , lower strain 
hardening and smaller rolling power requirement in general. 
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observed that the increase in densif i cation is more sensitive to 
change in friction factor at higher percentage reduction- This 
trend can be attributed to the larger arc of contact and higher 
relative density . The backward and forward slip are seem to 
increase with percentage reduction; in particular, the backward 
slip becomes very high in magnitude at higher reduction for all 
friction factors. As regards, the rolling power, both the 
friction and the plastic deformation components show an increase 
with P'er cent age reduction. While, friction work increases 
because of larger contact area, the plastic work increases due to 
more straining of material. 

From Fig. 5.9d it can be seen that the exit relative density 
increases linearly with percentage reduction for a given set of 
rolling conditions. However, at higher reduction (E5'/. and above) 
the exit relative density shows a decreasing tendency. This can 
be explained in terms of larger backward slip occurring at higher 
reductions. It is observed that the exponent (b> appearing in 
equation (S.44) decreases rapidly with increase in slip. This in 
turn, affects the densif ication and exit relative density falls 
at higher reduction. Also larger slip causes improper gripping 
of the strip between the rolls and thus leads to less 
densif ication. 

The numerical values of exit density are summarised in 
Tables 5.4 for different friction factors and percentage 

reductions. 

(iii) Roll Radius* The dependence of exit density and power 

shown in Figures 5.10a and 5.10b. 


components upon roll radius is 
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FIG. 5.10a EFFECT OF ROLL RADIUS ON DENSIFICATION 
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II IS evident that the densif ication behaviour is not 
significantly influenced by roll radius. This finding supports 
the conclusion derived from a microscopic view of a porous strip 
rolling, which indicates that the relative density is only a 
function of normalised contact area and slip between particles. 
Clearly both of these microscopic factors are not influenced by 
the distance travelled by the particles in the deformation 2 one. 
The exit relative density therefore is fairly insensitive to 
variation in roll radius. 

The rolling work, on the other hand increases slightly with 
roll radius. This is contributed partly by friction work and 
partly fay plastic work. The increase in friction work is due to 
larger arc of contact and that in the plastic work is because of 
the larger size of the plastic deformation zone. The results for 
variation of work and density with roll radius are given in 
Tables 5.5 and 5.6. 

<iv) Roll Speed: It can be from Figures 5.11a and 5.11b that 
the densif i cation is unaffected by roll speed, while power 
increases linearly with roll speed. The insensitivity of the 
densif i cation to roll speed is again related to the microscopic 
features of porous . metal rolling discussed in the previous 
section. The roll speed influences only the travel time of a 
particle through the roll gap and does not affect the 
interparticle contact area or slip between them. On the other 
hand, higher roll speed results in larger strain rates which in 
turn, leads to higher plastic deformation work. The backward 
slip also tends to increase with roll speed thereby causina 
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FIG 5.11b EFFECT-OF SOLt SPEED ON WORK 



larger frictional dissipation. Of these two contributions, the 
change in p'lastic def ormation power is more as compared to that 
of the friction power. These results are summarised in Tables 
5.7 and 5.8. 


(V) Strip Thickness: Figure 5.1Ea shows that larger 
densif ication occurs for smaller strip thickness. This can be 
attributed to the smoother entry of the strip between the rolls 
when it is thin. Consequently the entrance strip velocity is 
also larger. 

The plastic work increases steeply with initial thickness of 
the strip <Figure 5.12fa). This is obviously a consequence of 
larger material volume undergoing plastic deformation. The 
friction work also increases due to larger arc of contact and an 
increase in the magnitude of backward slip. The total rolling 
power, therefore, exhibits a sharp linear increase with initial 
slif' thickness. The results are shown in Tables 5.9 and 5.10. 

(vi) Initial Relative Density of the Strips The incremental 
change in relative density from the entrance to the exit sections 
is shown in figure 5.13a as a function of the initial relative 
density. It is observed that the percentage increase in density 
decreases for higher initial relative density. This is because 
of larger contact area and smaller available void space for a 
larger initial relative density. Due to these two factors, the 
resistance offered by the material for further deformation 

a hiflher yield stress value. On the 
increases resulting m a nigner yiHxu 

u .4 ti-a wnrk oer unit width increases with initial 
other hand the worK per utixv 
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re-lstive deiisity (Fig. 5.13b). Evidently the higher yield stress 
value increases the plastic work contribution. Also, because of 
larger actual contact area between the rolls and the porous 
strip at higher relative density as well as the higher yield 
stress value, the friction work will also increases. The results 
are shown in Tables 5.11 and 5.12. 


J5.4 COMPARISON BETWEEN THEORY AND EXPERIMENTS 


For multipass cold densif ication rolling, the exit density 
and linear dimensions of the strip were measured after each pass. 
A theoretical study of multipass rolling with conditions similar 
to those of the experiments was also conducted, in order to 
validate the present theoretical approach and the numerical 
solution procedure. The measured and predicted exit relative 
density values after each pass have been plotted against the 
cumulative percentage reduction in Fig. (5.14a>. It is observed 
that, there is reasonable agreement between theory and 
experiment, especially at higher cumulative percentage 
reductions. It is to be noted, however, that the actual friction 
factor (which has a significant effect upon densif ication) could 
not be measured and the value of 0.4 used for obtaining the 
results is only an assumed value. This value of friction factor 
seems to be a reasonable one for the present study in view of the 
fact that the roll srface was not smooth one and the rolling was 

performed under unlubriW, ted condition. Also, there ie a slight 
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is represented by the 
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From the measured dimensions of the strip before and after 
each pass the true length and thickness strains have been 
calculated. These are plotted against each other in Figure 
5.14fa. The figure also includes the corresponding theoretical 
predictions and the curve for a fully dense material undergoing 

ideal plane strain rolling. It is observed that the ratio ® w 

is less than one for porous strip, w'^ile it is equal to one for 
fully dense material. When the solid metal undergoes rolling 
de f o rmat i on under plane strain condition, the true strains in 

both longitudinal and thickness directions are same. However, 

when a porous metal is rolled under plane strain situation, most 
of the def oriTiat i on initially takes place in thickness direction 
which leads to closing of pores and thereby increasing the 
density of the strip. Very little elongation of the strip takes 
place. As the rolling continues further, a stage comes when the 
strip starts deforming both in thickness direction and 

longitudinal direction. When the porouus strip has been fully 
densified, then the true longitudinal strain ( Vs. tre 

thickness strain (Cy) curve becomes parallel to that of the solid 
metal. A good correspondence is observed between the 

theoretical and experimental curves shown in the figure. The 
experimental curve lies slightly above that of the theoretical 
curve. The following reasons can be attributed for this 
behaviour* (i) the existence of strain in the third direction 
( e,> i.e. width direction for the experimental situation as 
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FIG 5.14a COMPARISON BETWEEN EXPERIMENTAL AND 

theoretical RESllLTS FOR RELATIVE DENSITY 
FOR PWLWA5S ROLiING OF SINTERED COPPER 
STRIP 
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FIG.5.14b COMPARISON BETWEEN EXPERIMENTAL AND 
THEORETICAL RESULTS FOR STRAINS FOR 
MULTI PASS ROLLING OF COPPER STRIP 
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i.p fePtTi from Table 5.13 which leads to a reduction of the length 
strain <€■.,>; (ii) the actual friction factor under experimental 
conditions is not known and only an assumed value of 0.4 has been 
used in the theoretical cal culations . The results are summarised 
in Tables 5.13 and 5.14. 

5.5 ANALYSIS OF MULTI-PASS ROLLING WITH AND WITHOUT INTERMEDIATE 
ANNEALING: 

Figure (5.15), the path followed by the density curve 
as a function of reduction has been shown for multi-pass cold 
rolling with intermediate annealing. It has been assumed, while 
predicting the results shown in the figure, that relative density 
is unaltered by intermediate annealing. It is evident from the 
figure that, the incremental change in relative density for a 
given percentage reduction decreases after each pass. This is ■ 

because of the larger resistance to the deformation offered at 
higher material density, which was discussed earlier. 

A comparison between the exit density predictions for 
multi-pass rolling with and without intermediate annealing is 
shown in Figure 5.16. It is clear from the figure that the 
theoretically predicted values of the relative density of the : 

strip after cold densif ication rolling with and without j 

interiTiediate annealing are approximately same. fl similar ! 

observation has been found experimentally while hot rolling i 

sponge iron strip of density 1.9 git/c-c- at S50 C and 1130 C. ^ 

It was observed that the density is independent of hot rolling j 

temperature, while the ductility of the strip hot rolled at j 

o ». • t than that of hot rolled at 850°C OipO. A | 

1130°C was higher than that ot 

i 
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FIG 5-15 COLD DENSIFICATION ROLLING WITH INTERMEDIATE 
ANNEALING. ’ ' 
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siniilar trend was also found while cold rolling a copper strip 
With and without annealing C20D. From this present analysis it 
seems that the densif ication of porous copper strip is not 
affected by the mode of cold rolling , i.e. cold rolling with 
intermediate annealing or cold rolling without intermediate 
annealing. However, from the point of view of the properties 
achieved in the rolled strip the mode of rolling is an important 
f actor . 


A comparison of the incremental work for each pass for 
rolling with and without intermediate annealing is shown m 
Figure 5.17. ftUhough, densif ication is not significantly 
affected by annealing it is seen that »ork almost increases by 
factor of t»o, if annealing is not performed. This difference 
is because of work hardening effect. It is also observed that the 


worf p‘©f p*3ss 
process. This can be 
thickneisS as evident 


at later stages of the multi pass rolling 
attributed to the reduction in the strip 
from Fig. 5.1Eb. The results for cold 


iKnn ujith and without intermediate annealing are 
densif ication rolling with ana 

Tni-xat: s i*S and 5.16 respectively, 
shown in Tables 5.13 ana r 
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CHAPTER VI 
CONCLUSIONS 


6.1 CONCLUSIONS OF THE PRESfetfT STUDY 

(i) The eiaplified iipper bound technique based on streamline 
approach is shown to-be successful for the prediction of 
density and velocity fields and other parameters of 
interest in rolling. 

(ii) It is demonstrated that, there exists a global minimum 
value for rolling power which corresponds to proper 
dens i f i cat i on and metal- flow rate. 

(iii) The velocity, density and strain distributions indicate 
that compression is essentially homogeneous. 

(iv) The backward an’d ' forward slip and the roll pressure 
increase with percenttage reduction. The neutral point 
occurs closer to the esiit side than to the entry for all 
reductions. It is seen that the roll pressure is .such 
smaller as compared that of a solid material under similar 


rolling conditions. 

Exit relative density, increases with friction factor, 
especially at higher percentage reductions. The rolling 
power also increases with friction factor. With 

j v 4 i!i'nc-i tv reaches a maximuro with other 

percentage reduction density reacnes 

• inn the same. The power, on the other hand 
conditions remaining the same. i- 

j ♦ • atcn The roll speed and roll 
increases with reducti 

effect upon densif ication. The 
radius have insignificant ettecv k 

rolling power was found to increase linearly with roll 
speed. . Pdf d -igher initial strip density. the 



9l 

dfrisification is less anH mii- 

iess and rolling power is more. It also 

hp seen that thinner • 

thinner strips experience greater 

derisjf ication. 

<v.) Tt, coa,parison b.Ween theoretical predictions and the 
experimental measurements shows reasonable agreement for 
densification as well as strains. 

(VI 1 ) The densification trend remains same for multi-pass 
rolling with and without intermediate annealing. The work 
on the other hand is much smaller for annealed strip. 

SCOPE FOR FUTURE WORK 

^il ^ detailed experimental study can be conducted to verify 
the effects of various rolling parameters upon 
densification. 

(ii) The present results can be investigated by using more 

accurate FEM approach which takes into account realistic i 

friction conditions. 

(iii) The reason for the non-existence of global maximum of | 

I 

power consumed for higher percentage reduction in single | 

i 

pass rolling needs to be investigaed. i 

Civ) The study can be extended to hot densif ication rolling and 

also to include features such as front and back tension | 

, i 

in cold rolling. 

h ' 
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iBLE 5. 

1 : The 

Global ninimuuro 
Set of Rolling 

of Power 
Conditions 

Consumed for a 

given 




«! = 435 


«2 

= 445 



UF 

UP 

UT 

UF 

UP 

UT 

1 

3.415 

16.927 

21.126 

3.753 

17.316 

21.871 

.5 

3 . 322 

16.769 

20.87 

3.565 

17.155 

21.522 

6 

3.253 

16.625 

20.66 

3.396 

17.007 

21.20E 

.7 

3.219 

16.494 

20.498 

3.241 

16.874 

20.917 

.8 

3 .226 

16.377 

20.38 

3.110 

16.754 

20.66 

.9 

3 . 284 

16.274 

20.34 

2.999 

16.648 

20.44f 

.0 

3.413 

16.184 

20.382 

2.916 

16,556 

20.27E 

.1 

3.670 

16.109 

20.563 

2.869 

16.479 

20.15C 

.1 

4 .034 

16.047 

20.865 

2.859 

16.416 

20.07£ 

.3 

4 .400 

16.000 

21.183 

2.899 

16.368 

20.06? 


4 .767 

15. 966 

21.517 

3.006 

16.333 

20. 01^ 

.5 

5.136 

15.947 

21.867 

3.221 

16.314 

20.33' 

j 

.6 

5.507 

15.942 

22.233 

3.581 

16.308 

20.69: 

1 

E 

.7 

5 . 880 

15.951 

22.615 

3.943 

16.319 

21.06i 

i 

.8 

6.255 

15.974 

23.013 

4.307 

16 . 341 

21 .45C 

1 

.9 

6.633 

16.010 

23.428 

4.673 

16.378 

21.85^ 

1 

.0 

7.013 

16.061 

23.858 

5.041 

16.430 

22.27* 

.1 

7.395 

16.124 

24.304 

5.411 

16.495 

22.70| 

.2 

7.780 

16.202 

24.766 

5.783 

16.574 

23.16 j 

.3 

8.168 

16.292 

25.244 

6.158 

16.667 

23.62j 

.4 

8.558 

16,395 

25.737 

6.535 

. 16.772 

24.10 

'S 

8.951 

16.451 

26.247 

6.914 

16.891 

24.60 

*«» 

.6 

.T 

9.346 

16.639 

26.670 

7.296 

17.022 

25.12 

9.745 

16.780 

27.309 

7.681 

17.166 

25.64 



93 


2 . 8 

10.146 

16.932 

27.862 

8.067 

17.321 

26 . IS 

2.9 

10.550 

17.096 

28.430 

8.457 

17.489 

26. 7 ( 

3.0 

10.957 

17.271 

29.012 

8.853 

17.668 

27.31 








b 


Uj = 460 


'^2 

= 475 


UF 

UP 

UT 

UF 

UP 

UT 

0.4 

4 . 943 

17.899 

23.671 

6.739 

18.483 

26.07 

0.5 

4 . 680 

17.733 

23.242 

6.446 

18.311 

25.61 

0.6 

4.424 

17.581 

22.834 

6.152 

18.154 

25.16 

0.7 

4.179 

17.442 

22.451 

5.859 

18.011 

24.72 

0.8 

3.936 

17.318 

22.084 

5.571 

17.883 

24.31 

0.9 

3.705 

17.209 

21.743 

5.285 

17.770 

23.91 

1.0 

3.487 

17.114 

21.430 

5.005 

17.673 

23 . 53 ! 

1.1 

3.227 

17.035 

21.140 

4.727 

17.590 

23.17 

1.2 

3.083 

16.969 

20.882 

4.453 

17.523 

22.83 

1.3 

2.902 

16.919 

20.651 

4.184 

17.471 

22.51 

1.4 

2.743 

16.884 

20.456 

3.924 

17.435 

22.21 

1.5 

2.601 

16.864 

20.294 

3.668 

17.414 

21.93 

1.6 

2.487 

16.858 

20.174 

3.421 

17.408 

21.68 

I 

1.7 

2.404 

16.868 

20.101 

3.185 

17.418 

21 . 4 ^ 

1 

1.8 

2.355 

16.892 

20.076 

2.959 

17.443 

21.25 

1.9 

2.349 

16.930 

20.109 

2.747 

17.483 

21 . 0 € 

I 

2.0 

2.398 

16.985 

20.211 

2.548 

17.537 

20.94 

2.1 

2.523 

17.051 

20.404 

2.368 

17.607 

20 . 83 | 

2.2 

2.788 

17.133 

20-750 

J: . ' ■ 

2.210 

17.691 

20 . 7 e | 

2 . 3 

3.144 

17.228 

21.202 

i 

2.075 

17.790 

20 . 7 L 

2.4 

3.501 

17.338 

\ 'A ii, 

21.448 

1.970 , 

17.903 
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2.5 

3.860 


17.460 

22.149 

1.900 18.029 

20.78 

2.6 

4.221 


17.596 

22.646 

1.872 18.170 

20 . 8 S 

2.7 

4.584 


17.744 

23.158 

1.894 18.323 

21.07 

2.8 

4.949 


17.905 

23.684 

1.986 18.489 

21.33 

2.9 

5.317 


18.078 

24.225 

2.180 18.668 

21.70 

3.0 

5.686 


18.263 

24.779 

2.524 18.859 

22.22 
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Ul = 453 

* 




UF 


UP 

UT 


0.4 


4.297 


17-627 

22.741 


0-5 


4.064 


17.463 

22.344 


0.6 


3.844 


17.313 

21.974 


0.7 


3.630 


17.177 

21 .624 


0.8 


3.431 


17.055 

21.303 

( 

0.9 


3.243 


16.947 

21.007 


1.0 


3.075 


16.854 

20.745 


1.1 


2.923 


16.775 

20.515 

I 

1.2 


2.796 


16.711 

20.324 

1 

1.3 


2.690 


16.662 

20.169 

i 

1.4 


2.618 


16.627 

20.062 

1 

1.5 


2.581 


16.607 

20.005 

! 

1.55 


2.580 


16.603 

19.999 

1 

1.6 


2.290 


16.602 

20.008 

[ 

[ 

i 

1.7 


2.656 


16.611 

20.083 

[ 

1.8 


2.803 


16.635 

20.254 


1.9 


3.105 


16.673 

20.594 


2.0 


3,463 


16.725 

21.005 
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2.1 

3.824 

2.2 

4 .186 

2.3 

4.550 

2.4 

4.917 

2.5 

5.285 

2.6 

5.656 

2.7 

6.029 

2.8 

6.404 

2.9 

6.782 

3.0 

7.162 


16.792 

21 . 432 

16.872 

21.875 

16.966 

22.333 

17.074 

22.807 

17.194 

23.296 

17.328 

23.801 

17.474 

24.320 

17.633 

24.854 

17.803 

25.402 

17.985 

25.964 
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Tabl« 5.2 : 

Percentage 

Reduction 

1 

Effect of 

Friction Factor 

on Exit Relative Density 


Friction Factor (m) 


0.3 

0.4 

0.5 

0.6 


Exit Relative 

Density (s) 


5 

0.7241 

0.7250 

0.7259 

0.7267 

10 

0.7443 

0.7447 

0.7489 

0.7508 

15 

0.7564 

0.7635 

0.7682 

0.7737 

20 

0.7673 

0.7744 

0.7782 

0.7832 
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Table 5.3 : Effect of Friction Factor on Work 
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T^bl 0 5 * 4 : 

Effect of Percentage Reduction 
for Different Friction Factors 

on Exit 

Relative Density 

Percentage 

Reduction 


Friction 

Factor (m) 


C%) 

0.3 

0.4 

0.5 


0.6 



Exit Relative 

Density 

(P) 


5 

0.7241 

0.7250 

0.7259 


0.7267 

10 

0.7443 

0.7447 

0.7489 


0.7508 

15 

0.7564 

0.7635 

0.7682 


0.7737 

20 

0.7673 

0.7744 

0.7782 


0.7832 

25 

0.7677 

0.7763 

0.7797 


0.7825 

30 

0.7543 

0.7628 

0.7761 


0.7810 
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TABLE 5.5 : 

Effect of Roll 

Radius on Exit Reiative Density 

SI .No. 

Roll Radius 
(mm) 

Exit Relative Density 
(P) 

1 

67.50 

0.7744 

2 

80.00 

0.7744 

3 

90.00 

0.7744 

4 

100.00 

0.7744 

5 

150.00 

0.7817 

6 

200.00 

0.7892 

7 

350.00 

0.7999 


TABLE 5.6 : Effect of Roll Radius on Work 


Roll Radius 
(,mm) 

Uork/Uidth (N- 

-nun/s) 

UF 

UP 

UT 

67.50 

3.218 

16.800 

20.844 

80.00 . 

3.516 

16.742 

21.017 

90.00 

3.737 

16.707 

21.159 

100.00 

3.946 

16.678 

21.303 

150.00 

4.598 

16.796 

21.952 

200.00 

4.985 

17.051 

22.525 

350.00 

6.057 

17.494 

23.925 


7 



TABLE 5.7 : Effect 

of Roll Speed 

on Exit Relative Density 

SI .No. 

Roll Speed 

C fflffl / e ) 

Exit Relative Density 

(P) 

1 

400.00 

0.7720 

2 

508.68 

0.7744 

3 

525.00 

0.77126 

4 

550.00 

0.77126 

5 

650.00 

0.77126 

6 

800.00 

0.77126 


TABLE 5.8 .* Effect of Roll Speed on Work 


SI .No. 


1 

2 

3 

4 

5 


Roll Speed Uork/Uidth (N-aim/s) 

(mni/e) 



UF 

UP 

UT 

400.00 

2.581 

13.161 

16.388 

508.68 

3.218 

16.800 

20.844 

525.00 

3.437 

17.228 

21.512 

550.00 

3.576 

18.071 

22.535 

650.00 

4.211 

21.370 

26.632 

800.00 

5.205 

26.282 

32.779 


6 



Tabl* 5.9 ; 


Initial Strip Thickness on Exit Relative 
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SI .No. 

Initial Strip Thickness 
(fflin) 

Exit 

Relative Density ( P) 

1 

5.0 


0.7406 


2 

4.0 


0.7406 


3 

3.0 


0.7419 


4 

2.0 


0.7447 


5 

0.5 


0.7571 


Table 

5.10 : Effect of Initial Strip Thickness on Uork 


SI .-No. 

Initial Strip Thickness 
(otin) 


Uork/Uidth 

CN-nun/a) 




HP 

UF 

UT 

1 

5.0 

17.831 

1.095 

21.380 

2 

4.0 

14.008 

0.994 

16.758 

3 

3.0 

10.327 

0.837 

12.306 

4 

2.0 

6.761 

0.636 

8.038 

5 

0.5 

1.706 

0.214 

2.000 



lU^ 


Table 5.11 


Effect of Initial 
Densif ication 


Relative 


Density 


of 


Strip 


on 


SI -No. 

Initial Relative 
Density of Strip 
( 

Exit Relative 
Density ( P) 

Percentage 
Increase In 
Relative 

Density (%) 

1 

0.70 

0.7744 

10.630 

2 

0.80 

0.8566 

7.084 

3 

0.90 

0.9299 

3.322 

Tatblife 5* 

12 : Effect of Initial Relative Density of 

Strip on Uork 

SI. No. 

Initial Relative 
Density of 

Strip (Pi) 

Uork/Uidth 

(N-mm/s) 


UP UF 

UT 

1 

0.70 

16.800 3.218 

20.844 

2 

0.80 

24.861 5.725 

31.706 

3 

0.90 

35.122 9.526 

46.102 








. y;: ^ 

% 

% o & ^ 


^ 



TABLE S.13 ! Experimental Results for Hulti-pass Rollinfi of Sintered Copper Strip 
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T*bl« 5.14 


Thsoretical Results 
Copper Strip 


for Multi-pass Rollinfi 


Sintered 


Nu 0 iber 
o ( the 
Pass 


Percentage Relative 

Thickness Density 

Reduction (p ) 

with Respect 
to Original 
Thi ckness 
<%) 


Length 
Strain 
(fix ) 


Thickness 

Strain 

(fiy ) 


Width 

Strain 


(Ez) 


1 

17.09 

0.8095 

0.12613 

-0.18741 

0 

2 

24.6 

0.8393 

0.1854 

-0.2826 

0 

3 

29.116 

0.8506 

0.2335 

-0.3440 

0 

4 

34.5 

0.8748 

0.2852 

-0.4238 

‘ 0 

5 

45.2 

0.9079 

0.4256 

-0.6013 

0 

6 

54.83 

0.9318 

0.5837 

-0.7855 

0 

7 

61.4 

0.9509 

0.7316 

-0.9536 

0 


I 


! 




15; Multi-pass Cold Densification Rollinfi with Intermediate 
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Y1 -STRIP THICKNESS AT THE ENTRY SECTION 
UO-STRIf VELOCITY AT THE ENTRY 


U-ROLL SPEED 
FR-FRICTION FACTOR 

S-EQUIVALENT STRESS FOR THE SIVEN MATERIAL 

YEOLD-STRIF THICKNESS AT THE EXIT 

EPS-C0NVER6ENCE CRITERION 

Ny-NUHBER OF ITERATIONS 

EL-UENTH OF CONTACT 

UORKI-TOTAL WORK FUNCTION 

RHOO-DENSITV AT THE ENTRY 

X81-STREAM LINE COORDINATE 

R-RA0IUS OF THE CIRCLE 

RA-RADIU8 OF THE ROLL 

RHO-DEN5ITV IN THE DEFORMATION ZONE 

U-AXIAL VELOCITY 

V-TRANSVER8E VELOCITY 


X-X-CORDIHATE 
Y-Y-CORDINATE 
8 I -HASS FLOW RATE 
XNODE-X-CORDINATE OF THE NODE 
YNODE-Y-COR0INATE OF THE NOOE 

RHONO0-DENSITY AT THE NOBE , 

UNOOE-AXIAL VELOCITY OF THE NODE 

VN00E-TRAH8VER8E VELOCITY OF THE NODE rn»nTWATE 

ui-OFRATIVE OF THE STREAM FUNCTION WITH RESPECT TO Y-CORDINATE 
W THE STrIaM FUNCTION WITH RESPECT TO X-CORDINARE 
FRICT-UORK loss due TO FRICTION 
DI8CT-W0RK LOSS AT THE VELOCITY DI8C0HTINUTY 
DEFORM-WORK REQUIRED FOR PLASTIC DEFORMATION 
FF— DERATIVE OF WORK FUNCTION WITH RESPECT TO Y20LD 
DFfSsECOND DERATIVE OF WORK FUNCTION WITH RESPECT TO YEOLD 


IMPLICIT REAL • 8<A-H,0-Z» 


RA-fi7.5 

Y1-I .0 

Y8-6.80 

FR-0.4 

UO-450 . 00 

UX-508.68 

RHOO-0.70 


EPS-0. OOOt 

BHAX-S.O 

DELB-O.OOO! 

DB-0 . I 

Bl-O.I 

HUH-50 


CALL «ORKtRHOO,Yl ,U0, Y£.FR,8A,UXxBl .WORK* > 
WORKIl-WORKI 

Bit— Bl*BEiiB . , Bi I uoHKl) 

CALL WORK ( RHOO , Y I , UO , Y8 , FR . RA . UX ' ® ” • 

UORKIIt-MORKl 





30 

0243 

to 

00 

C 

20 

2450 

C 


CAtt t^ORKCRHOO# Y1 .UO^YS FR RA 

woRKiie-yoRKi 


UX,B1S,yORKl ) 


FF 1 » « WORK 1 1 ! -WORK 11 2 )/ ( £ . 0#DELB) 

OF F > - ( WORK I n +W0RK 1 1 8- < £ . 0*UORK 1 1 ) ) / < DEUB. *2 ) 


B8*Bt4DB 

B2J-BB+DELB 

WORK 1 21 - WORK I ^ ” 
B88-B8-0EtB 

CALL WORK(RHOO. Y1 .UO, Y8,FR,RA,UX,BS8,U0RKI) 
MORKISE-WORKI 

FF8-{WORKI81-HORKI88)/<8.C*DELB) 
DFF8«(W0RKI814W0RKI8S-<8. 0*W0RK 1 8 ) > / ( BELB* *S ) 
IF (FF1*FF8) 10,80,30 


Bt-B8 

FFI-FF8 

IF<B8.6T.B«AX)S0 TO 9843 
60 TO 689 
TYPE », 'QUEST ' 

60 TO 84S6 


PAUSE 'IN TO ITERATION* 

BOLD-Bl 

FF-FF1 

OFF-OFF 1 

ITERATION-1 

BNEU-BOLO~ ( FF/DFF J 

CALL WORK tRH00,Yt,U0,Y8,FR,RA,UX, BNEU , WORK 1 1 

UORKIN-UORKI 

BX1-BNEW4DELB 

CALL WORK ( RHOO , Y 1 , UO , YS , FR , R A , OX , BX 1 , WORK I ) 

WORK IN- WORK I 
BXS-BNEW-OELS 

CALL WORKIRHOO,Y1,UO,Y8,FR,RA,OX,BX8,WORK1) 

UORKIt-UORKI 

FF- { WORK I H-UORK I L 1 / ( S . 0*DELB ) 
OPF-(UORKlN4UORKIL-t&.0*UORKIN)>/(OELB*«8) 

WR I TE { S , • > ' I TERAT I ON- ' , ITERAT I ON 

IF { DABS t FF J . LE . EPS ) STOP 

IFCBFF.LT.0.018TOP 

ITERAT I ON- 1 TERAT ION+ 1 

IF< ITERATION. 6E.NUH)ST0P 

BOLO-BNEW 

60 TO 60 

BNEW-BE 

STOP 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SOBROUT I HE WORK ( RHOO , Y 1 , «0 , re , FH , R A , OX , B . WORK n 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXWCXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE TO CALCULATE KNSITY, VELOCITY FIELD IN THE DEFORMATION 
ZONE AND WORK REQUIRED FOR FLA5T1C DEFORMATION FOR STRIP ROLLINQ 
UNDER plane strain 


EXTERNAL SUBROUTINE 

SHAPE# 



C - 04ACOB 

C “ DRIVES 

C - RLRStIC 

C 


IHPLICIT REAL • 8(A-H,0-2) 


DIMENSION P0SeB(E},UEieB(&),SHAPI(4) 

I ,x<n,EE),ytii,ss),x5i(n,eE),si(n,EE) 
e ,R(it,SE),Ftii,eE>,z<it,ee),RH 0 (it,eE>,u(ii,e 2 ),v(ii,eE> 

3 .RHOIteS) ,YI(BZ),XI(ZB),RI(2e) ,EQS( 1 1 , BE) , AVSY1 11 . EE) 

4 ,SXJ810,4),STtE10,4),SXYCS10,41,EQSRt11 ,eS),AVSXin,SS) 

6 ,EX(1 1 ,EE3 ,EYt n .EEl.EXYtll ,EE),EQSS(EtO,4),AV8XT(11,EE) 

7 ,ROLLn,SE),VELO(n.EE) 

COMHON/FEM/XNODE (210,4), YNODE (210,4) 

COnnON/CEH/UNODE ( S 1 0 , 4 ) , VNODE ( E 1 0 , 4 ) , RHONOD ( El 0 , 4 ) 


C 

c 


c 

c 

c 

c 


908 

C 

C 

c 


M 18 NO. OF CRID POINTS IN X-DIRECTION 
N IS NO. OF ORID POINTS IN Y-DIRECTION 
N»1 1 

n»EE 

EL-DSORTI (RA-({Y1-Y2)/2.0))*2.0»(Y1-YE) ) 
DXSI"1 . 0/DFLOAT(N-1 ) 

DX-EL/OFUOAT(M-1 ) 


EVALUATION OF DENSITY AT THE SURFACE OSINS DENSITY 
EXPRESSION 


TXVAL»-DX 
00 908 I"1,E8 
TXVAL-TXVAL+DX 
XI(I)«TXVAL 

RKD-I (EL««S)/(£,0*(Y1-YE)>) + ((Y1-Y2)/B.0) 
YR-R1(I)7YB 

YI ( 1 > —DSQRT t R 1 { I ) ••£- ( XI ( 1 ) -EL ) •*£ ) +YR 

RNIT-( 1 .0-RHOO)/RHOO 

VARl-«YltI)/Yt)»4B 

RHOKI )-1 .0/(1 .0+(RN1T*VARI)) 

IF(RH0I<I).6T.0.998)RH01I1)-C.998 

CONTINUE 

evaluation of density in the DOMAIN 


X8IVAL—DX81 
DO 100 .N 
XVAL*-DX 

XSIVAL-XSIVAL+DXSI 


0 ISO K«1,M 
VAL*XVAL+DX 
:<4.K)-XVAL 
;SI(J,K)»X8IVAL 

ufxiSi (£*•*£)/(£. 0*X8I(J,K)*(yi-Y2)))*(X81(J.K) 

-(Y1-Ye»/E.O) 

0-R<J,K>+XSI(J,K)*Ye 

DS<IRT(R<J.K)**2-tXtJ.K)-fL)**2)+Y0 

( I ( J . K 1 -RHOOfUO* Y t •X8 1 ( J . K ) 

1»Y(J.K)-Y0 

■2»Y2*{Y(4,X)-Yes 1 

•S-{Yt-YE)/£,f ,, , . 

r I -tl / ( TB+ { T|» (TJI+TS? ), ) 

■6-X(J.K)-EL 



VE«T4/<TS3:*KT3#(T4+T5>) ) 

F< J,K)«RHOO»UO*Y!*V1 
2 ( J , K ) •“RHOO*UO* Y 1 ♦Vfi 
RHOCJ,K)«RMOHK) 
il( J,K)»F<jr,K)/RHO<J,K) 

V( J,K)«Z{J,K>/RHO<J,K) 

CO TO 33 

42 VtJf,K)-0 0000 

RHO<J,K>«RHOI(K) 

33 8I<3,K)*RHO{ J,K)‘>U<J,K>*Y<J,K) 

ISO CONTINUE 

too CONTINUE 

J-t 

DO 358 1*1,88 
358 CONTINUE 


530 

663 

C 

C 

c 


1800 

1830 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


DO 663 

DO 530 K*1,£8 

VELO< J,K>*DS0RT(UCJ,K)**8+V(J,K)*»£) 

CONTINUE 

CONTINUE 


HUHBERINC THE NODES AND SETTING THE CORD 1 NATES 


NEL»0 

DO 1838 J»1,10 
DO 1800 K«1 ,81 
NEt*NEL«ft 

XNODC<NEt,1>»X(J,KJ 

YNODE(NEL,1)«Y(J,K) 

RHONODtNEtiT )«RHO(J,K) 

UNOOECNEt,tt*UCJr,K> 

VN0DE<NEL,1)«V(J,K) 

XNODE C NEL , E > «X ( J , K+ U 
YNODE<NEt,8)«Y<4f,K+t ) 

IIH0N0D(NCL,8)«RH0U#K+1 > 
yNODE(NEL,S)»U(J#K+l J 
VNODECNeL,C1»VC Jf,K+t ) 

XN0D8tNEI-,3>»XIJ+1 tK+l ) 

YH0DE<NEL,3>*YCJ-»*t#K+1 > 
RH0N0D(NEt#3I»*8H0lJ+1/K+1 ) 

UN0DE(NEt,3)*U(Jr+t »K+t > 

VN0DE<NEt,3>»V( J+1 iK+l ) 

XN0D£<NEt,4)«XCJ+1 ,KI 
yNODE<NEt,43*Y<J^1 ,K) 

RHONOD C NEL , 4 ) »RHO < 1 . K > 

yN0DEINEt,4)*U(J+1#K) 

VNODECNEL,43*V(J+l,K) 

CONTINUE 

CONTINUE 

WORK FUNCTION FOR PLASTIC DEFORHATION 


« STRAIN IN X-DIRECTION 

« strain RATE IN X-DIRECTION 

SX« AVERAGE STRAIN RATE IN X-DIRECTION 

• STRAIN IN Y-DIRECTION 

• STRAIN RATE IN niRECTXON 

SY* AVERAGE STRAIN RATE IN Y-DIRECTION 

:Y» shear STRAIN 

:Y» SHEAR STRAIN RATE 

^SXY» AVERAGE SHEAR STRAIN 

fSR* EOUIVALENT STRAIN 

(«« equivalent strain 


t 




c 

c 

c 


in 

c 

c 

c 


T03 

T02 

C 

C 

C 


704 

C 

c 

c 


70S 

C 

C 

C 


70S 

C 

C 

C 


707 

C 

C 

C 


EVALUATION OF STRAIN RATES IN THE DOMAIN U8INC STRAIN SUBROUTINE 

00 in iELEn»nEio 

NELCM»EtO 

NN0DL«4 

N00E«4 

CALL STRAINCNELEM, IELen,N0DE,NN0DL,SX,SY,SXY) , 

CONTINUE 

EVALUATION OF AVERACE STRAIN RATE AT A NODE INSIDE THE DOMAIN 
DO TOE 

DO 703 K»nEO 
1*1 J-t >* 2 nK 

AVSX( J+1 /K+l J»<SX< I,3I+5X(I+1 ,4)+SX( I+E1 ,E)+SX( I+EE;! ))/4.0 
AVSY< J^nK-n )»(SY(I,3>+SY<I+1 ,4l*»'SY(I+2t ,E)+SY<X+2E, t ))/4.0 
AVSXYI J+nK+1 )*(SXY( I,3)+SXY(I-H ,4)+SXYCI+2l ,£)+SXY( I4£E, 1 J ) 
/4.0 

CONTINUE 

CONTINUE 

EVALUATION OF AVERAGE STRAIN AT A NODE AT THE CENTRE LINE 


DO 704 K*neo 

AVSXU,KTt )-(8X(K,£)+SX<Kfin >)/2.0 
AV$Y( J,KT| >«{SYCK,£)+SYCK+in )>/e.O 
AVSXYC J,K+1 )*<8XY(K,£)+SXY<K+1 , 1 ))/2.0 
CONTINUE 

EVALUATION OF STRAIN RATE AT THE SURFACE 


J *10 

DO ros K^neo 

!•( j-t 

AVSXU 41 iK 4 n-< 8 X{I, 3 )+SX(I+ 1 , 4 n /£.0 
AV 8 Y<J’*‘t,K 4 l l«tSY<I, 3 )+SY{I+l, 4 )>/ 2.0 
AVSXYU+I ,K4! J-<SXY(I,3n-SXYn+l ,4))/2.0 
CONTINUE 

EVALUATION OF THE STRAIN RATE AT THE ENTRY SECTION 


<«t 

70S j-ns 

iv$x<3+t ,xi«<sx< i,4)4Sx<i+2i , n)/2- 0 

^VSYCJ4t ,KI*(SY(I,4)+SY(I+2tn >>/2.0 

4 VSX Y 1 3 + n K > » t SX Y 1 1 , 4 ) +SXY < I +£ t , D > • 0 

CONTINUE 

EVALUATION OF STRAIN RATE AT EXIT SECTION 

K*et 
DO 707 

AVSX< J4| |K+I 

AVS Y < 1 . K4 n • C S Y { 1 , 3 ) r £ ) ) /£ 0 

AV$XY(3+t ,K+t )*(SXY(I#3>+SXY(I+2n2>>'2,0 

CONTINUE 


:VA«-UMION ®F strain RATES AT THE CORNERS OF THE DOHAIN 


'SX< 1 »t>*‘SX{ 1.1 > 

SYff.IJ-SYd,?) 

< 



c 

c 

c 


S33 

6t1 

C 

c 

c 


615 


5B1 

416 

353 

C 

C 

C 


8600 

8638 

C 

C 

c 


4VSXCII , n»sx<i$o,4) 

4VSYC11 *! )«SY<t90»4) 

AVSXY<n,1>*8Xy(190,4) 

AVSXCI,86>*SX(81,2) 

AVSVCt .66J»SY(2t ,6) 

AVSXYCt ,68>«SXY(61,6) 

4VSXCttt .66)»SXC610,3) 

4VSY<tt,efi)«SY<S10,3) 

4VSXYC H ,62>«8XY(610,3) 

EVALUATION OF EQUIVALENT STRAIN AT NODES 

00 Sit 3-1,11 
00 633 K»l,66 

TT1««AVSX( J,K)-AVSY(J,K))**2 
TT2«AVSX<J,K)+AVSY(3,K) 

TT3-(AVSX<J,K))**2 
’rT4-(AVSY(J,K))**2 
TTS-<RH01 .5 

TT6-1 .0/<2.49*(l ,0“RHO(3,K))**0.5U) 

TT7«(AVSXY{J,K) J»*2 

EQSR ( J , K ) -TTB^DSQRT t U 2 . 0/9 . 0 ) * < TT 1 +TT3+TT4+ 1 :|*^*TT7 ) ) ) 

CONTINUE 

CONTINUE 

EVALUATION OF CUHULATIVE STRAINS 


00 3S3 

00 416 K-I,E2 

IF<K.EQ.1 ISO TO SIS 

AVEQSR - 1 EQSR U , K- 1 1 +£086 ( J , K ) ) /2 . 0 

AVEX- < A VSX ( J , K-n + AV8X < Jf , K ) ) /6 , 0 

AVE Y* ( AVS Y < jr , K- 1 > + AVS Y ( J , K ) ) /£ . 0 

H VEXY- < A VSXY ( 3 , K- n ^AVSXY ( 3 , K ) ) /6 , 0 

AVUVAL«(UC3,X-'1 )'*'UC3,K) 1/6.0 

lQS(3,K>»EQ$(3fK-! I*KAVEQ8R*( <X( J,K)*X<3,K-t ) J/AVUVAL) ) 

£XC3,K>»£X<3,K-'1 >^(AVEX*C (X(3,K1-X(3,K-1 ) 1/AVUVAL) > 

EY<3#K1-EY13#K“"1 )+<AVEY*< <Xt3|K)-Xt3,K-l) )/AVUVAL) ) 

CXYt3,K)-EXYI J,K-1 HCAVEXY*! (X(3,K)-X(3,K-1 ) )/AVUJ^AL) ) 

80 TO 561 

EX<3,Kl-0.0 

EY(3,K>«0.© 

E88(3,K>-0.0 

EXYC3,IO«0.0 

CONTINUE 

CONTINUE 

CONTINUE 

ttUMBERINC THE NODES FOR EQUIVALENT STRAIN 


lELEH-0 

00 8638 3-1,10 

©0 8800 K»t,6t 

IELEH-IELEH+1 

EQSS i 1 ELEH , 1 ) »EQS< 3 , KI 

EQ8S< IELEH,6)»E08^3,X+I J 

£QSS< I ELEH , 3 ) -EOS ( 3+ 1 , K+U 

£8SS{ IELEH,41-EQS{ 3+1 ,K> 

CWTINUE 

CONTINUE 

EVALUATION OF ROLL PRESSURE 


3-1 1 

00 8309 K-t,66 



C 

C 

c 


641 

C 

c 

c 


c 

c 

c 


8 $a 

c 

c 

c 


PI 1*431 . 004UEO$<J,K1 + O.O1 )*»0.32) 

PIS*Ct O^AVSrU.K) )/3.0 

ro-l .0-PHO{J.KJ)**0.514) 

PI4*8.O/C0.0*FO**8) 

pi 5 *ro*»« 

P16»6VSX( JiKI-f-AySrCJ^K) 

PI7*<P10«P! 1 )/E08RU,K) 
P1d*(E.0*AV$X(J,K))/3.O 
PI$«Pl7<.CPie+< ( 1 .0-P14)*P!5*P16)) 
P«0*r!74<Pie^< < 1 .0-Pt4)*PI5*Pl6)) 
POLL < J , K > "DSCIPT < P 1 84*e+PeO*4S ) 
COKTIKUE 


EVALUATION OF PLASTIC WORK 


F|«*0 . 0 

00 641 IELE«*1,£t0 

NELE«*2I0 

N«O0L*4 

NOAUS-E 

CALL 0BIVE8<NELErf,EQSS, 

I NCAUS , NNO0L , SUH , lELEH ) 

FI>*FI+8UH 

CONTINUE 

OEFOPH*FI 

WORK LOSSES AT THE ROLL SURFACE AND AT THE DISCONTINUITY 


F066B< I J-0.5T7350E69S 
P0S6B(2)»-P0S6B{1 ) 

UEXSB( 1 )"1 • 0000000 
HElCB<e»-1 .0000000 
SHAPKI »■< 1 .0-POS6B< 1 ) )/2. 0 

8HAPI<B»-n .0+POS6B(n)/8.0 
SHAPIOJ"!! . 0-P086B<fi> )/B. 0 
SHAPIt4>“<J .0+P086B<E))/8.0 

WORK FUNCTION AT TWE ROLL SURFACE 


CO . 0 
J-I1 

DO 888 K»i.8t 

DDX»XlJ,K-H»-XtJ,K) 

DDY*Y<aiKJ-Y(<I.K+l ) 

DDL-D8QRT { DDX**2+I>DY**8 » 

EcfsHAPI < 1 )*RH 0 ( J>K)+SHAPI ( 8 )*RH 0 ( J»K +1 * 

E2»(SHAPI C3)*RHO<3»K)+8HAPI tAJ*RHO< J«K*1 > ***^'^o> 



i»r|s 43 l O 04 ((EFt'«'O.OI}**O. 3 E) 

C!»£|40Afi8<UX-01 K+1 ) 

EFfi-eHAPl < 3 J -EOS t a , K )+SHAPn 4 >*EOS< J , K+ 1 ) 

RF2»43t .00*ttEF8+O.On**6.32) 

CS>'E8*0AB$ < OX-BE ) oWElOB 1 8 )*DET JL*RF8 

CONTINUE , 

FRlCT«(6*F8J/tSQRT<3.0)) 

UORK FUNCTION AT THE DISCONTINUITY 


0 * 0*1 





444 


C 

C 

C 

C 

c 


c 


K «1 


00 444 . 10 

DFY»Yt Jf + I ,KJ-Y<J,K) 

DETJH-DFY/fi.O 

E3-<8HAPHn*RH0t J,K)+8H4PI(8)*RH0(>I+1 ,K))»*8.S 
E4>(SHAP!(3)*RH0( J,K)')’8HAPI(4)*RH0(J+I ,K))**8.5 
03«B80RT<«8HAPI(1 )*V< J,K)+8HAPI{8J*VtJ+l .K>)*«8) 

04»D8QRT< (SHAPI(3>*V<3,KJ+SHAPI<4)*V<J+1 .K))**S) 

EFS-SHAP 1 Cl > ‘EOS U , K ) ♦SHAP 1 1 8 ) *EOS C J+ 1 , K ) 

RF3-431 .0©*C(EF3+0.01)**0.38) 

C3oE3»03*UElCB(t >*DETJH*RF3 
EF4-8HAPH3)*EQ8(J,K)+8HAPI(4>*E0S(J+1,K) 

RF4-431 .00*((EF4+0.01)**0.38) 

C4«E4*D4*UEICB(8»*DETJH*RF4 

Q«Q+tft3+C4» 

CONTINUE 

BI8CT-8/{80RT{3.0)> 

UORK I -DEFORW+FR 1 CT+DISCT 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


SUBROUTINE SHAPE4tDERIV,SHAPE,XE(JlV, YEQIV) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DERIV-OERATIVE OF SHAPE FUNCTION UlTH RESPECT TO LOCAL COORDINATES 


XEQIV, YEOIV-LOCAL COORDINATES 
X.Y-LOCAL NAME FOR XEOIV AND YEQIV 


IMPLICIT REAL * 8tA-H,0-Z) 
DIMENSION DERIVCE.ei.SHAPECS) 


X-XEQIV 

Y-YEQIV 

XY-XPY 

SHAPE C 1 J - C 1 . -X-Y+XY J * . 8S 
SHAPE t E >- C 1 . +X-Y-XY J * . 85 
SHAPE ( 3 ) - t 1 . ♦X+Y+XY ) * . 85 
SHAPE(4J-<1 .-X+Y-XY)*.85 
DERIVCI.l J-C-1 .♦YIP. 85 
0ERIV«1,8>«C1 .-Y>*.85 
0ERIV<1.3)i-(t .♦TJP.ES 
DERIVCI .4)-t-l .-Y)P.85 
DERIV(E.U*(-1-+X)*.85 
0ERIVt8,8tAC-1 .-XJ4 .es 
DERIV f8.3>"Cl .+X)4.85 
DERIVCE,4J«C1 .-XJ4.8S 


RETURN 

SSxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE D JACOB t DERI V . DET JB , D JACK , D JAC I , NNODL 

iJSxxixSxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DETJB-DETERIMINANT OF JACOBIAN MATRIA 
DJACI-INVERSE of JACOBIAN MATRIX 
0 JACK- JACOBI AN MATRIX r.EMENT 

NNODL-NUMBER of NOBEf i" 

Mei.£H«iNUHB£R OF ELEHENT8 IN WESH 


*^N8 1 QN**DER I V C 8 , 8 ) . D J ACI 1 8 . 8 > , D J ACK < 8 . 8 ) 
dSoH/FEM/XNOBE < 81 0 , 4 > . VNODE (810,4) 

T ,» tPmeoBARY MATRIX TO ALLOW THE JACOBIAN 


TO BE FORMED 


TEWrt!«" 0.0 
TEffFYE»0 . 0 
TEPIFT3-0. 0 



so 


c 

c 

c 

30 


40 

C 

C 

C 


C 


C 

C 


1 

e 


c 

c 

c 

c 

c 


4444 

4445 

4446 

4447 
48 &S 


TErtPY4»0 0 

00 SO IN0DL»1,4 

XNODI-XNOOEf lELEH, INODL) 

yi«ODI*YNO0€< lELEfl, INODLJ 

TEHP Y t •TEHP 1 1 +0ER 1 V (1 , INODL ) *XNOD I 

TE«FyE»TE«PrE*»*0ERIV( T , IWOOLI*yNODI 

TEriFY3«TE«Py3+DERlV(e,IN0DL>*XN00I 

TEi1Py4*TEHFY4+DERIV(£*IN0DL>*YN0DI 

CONTINUE 

0JACKC 1 , I )»«TEHFY1 
DJ4CK( T ,S)»TEFFy2 
0JACK<e,1 >«TEHPY3 
0JACK<E,2)*«TE«PY4 

DETJB«*0JACK( 1 , t >»DJACKt2,S)-0JACK( 1 ,e)*DJACK{£> 1 ) 
CHECK FOR ZERO OR NECATIVE DETERMINANT 


IF<0ETJB)3O, 30,40 
CONTINUE 

WRITEtS^eOOOUELEM 

STOP 

CONTINUE 

INVERT TEMPORARY MATRIX TO FORM JACOBIAN 


OJACI (1,1 )*0JACK(E,E)/0ETJB 
0 JAC I « E , E ) *0 JACK (1,1) /DET JB 
0 J AC 1(1,2 ) «-0 J ACK (1,2) /DET JB 
0JACI (2, 1 )»-0JACK(2, t )/0ETJB 
FORMATOX, 'IELEM*',14) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE STRAINtNELEH, 1ELE«,N00E,NN0DL,SX,SY,SXY» 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


IMPLICIT REAL * 8(A-H,0-2) 

DIMENSIOK DERIV(E,8J.8HAPE<8J, 

0JACUE,e>.nJACKte,E).CARTL(a,4J, 

8XtEt0.4),#Y«810,4J.8XY<C10,4) 

COMMON/PEM/XMOOE 1 8 1 9 , 4 ) , YM0DE< 81 0,4) 
COMMON/eEH/UNODE (810,4), VNODE (810,4), 


RHONOD(810,4) 


8X<IELEM,NODE)“ 

SYdELEM.MODE)- 

8XY<IELEM,N0DE> 


X-8TRAIN RATE AT NODE 
Y-STRAIN RATE AT NODE 
■ SHEAR-STRAIN RATE AT NODE 


10 545 NODE- 1, 4 

;f{MO0E.eq.i ) CO to 4444 

;F(N0DE.EQ.8> 60 TO 4445 

[FtNOOE.EO.S) CO TO 4446 

[F(NO0E.EO.4) CO TO 4447 

CEQIV— 1 .0 

rEQIV— 1 .0 

)0 TO 4S88 

(EftlV-l .0 

fEQIV— 1 .0 

M) TO 4888 

<EQ1V-1 . 0 

YEOIV-1 . 0 

BO TO 4888 

KEOIV— 1 .0 

YECiIV"1 .0 

BO TO 4888 

S?I‘J,*^!Br^ti»PRIV.SHAPE,XEOIV.YEQIV) 



C^LL DJACOB(DERIV,DETJB,DJACK,DJACX iNNODL, lELEWiNELEK) 

DO 300 I0l«E-1,e 
DO 300 INDDL«t,^ 

CARTt<IDIHE,INODL)»0.0 
DO $34 JDIHE»I,8 

CAliTLC!DIH£,lNODL)*CARTtCIDl«E,lNODL)+DJACI(IDI«E,JDIHE)* 

1 DERIV(J0I»E,INODL) 

934 CONTINUE 

300 CONTINUE 

DUD)(«0 . 0 
DUDY-0.0 
DUDX-0 . 0 
0VDVO.O 

DO 489 lN0Dl.-t,4 

DUDX-DUDX+CARTLt t , 1N0DL)*UN0DE< lELEH, INODL) 
DUDy*DUDy4.CARTtC8,IN0DL>*UN0DE(IEtE«,XN0DL) 

DVDX*DVDX-i-CARTL< t , INODt )*VHODE( lELEH, INODL) 
0VDY*DVDY-»-CARTL<8,INODt)*VNODE<IELEM, INODL) 

489 CONTINUE 

SX(IELEH,NODE)-DUDX 
SY«IELEH,NODE)*DVDY 
$X Y < I ELEH , NODE ) •DUDY+DVDX 
S45 CONTINUE 

RETURN 
END 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE 0RXVE8(NELEH,EQSS, 

I NCAUS , NNODL , SUN , lELEW) 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c W 0 AUS*NU«BER of CAUSS POINTS IN RULE 

c POSCLCNCAUS) -POSITION OF GAUSS POINT IN LOCAL COORDINATES 

C «EICLtNGAUB>-CAUS8 WEIGHING FACTORS AT SAMPLING POINTS 

C AREAW-PRODUCT of ELEMENT AREA AND GAUSS WEIGHING FUNCTIONS 


C 

C 

C 

C 

C 

c 


c 

c 

c 


c 

c 

c 


external subroutine 

SHAPE4 

OJACOB 

PLASTIC 


IMPLICIT REAL ♦ S{A-H,0-Z) 

DIMENSION DERIV<8iS>*AREAW<4),SHAPEt8), 

DJACnSiB) .D3ACKtBiB)#CARTL<B,4), 

POSGLtB^ 

WEICL<8)iW(4),EGSSI2f0.4),8UMY<4), 

SUMt C 4 ) , SUME C 4 ) , SUM3( 4 ) , SUM4(4 ) , SUMI 1 4 ) 

COHHOH/FE»/XHOBE<eiO,4),YNODEt210.4) 

C0WM0N/6EH/tlM0DE {S10.4),VNODE(81O,4>,RHOHODt£10,4> 

SET UP P08IT10M8 AND UE1CHT8 FOR 2 POINT 8AU8S RULE 

POS6L(t )-0. 5773502 
P086L(E)— POSCLCt) 

IffilCLO >-1 .0 
»ffilSLte)-WEISL<1) 

CALCULATE SHAPE FUNCTION AND DERIVATIVES FOR ELEHENT8 


LCAU8-0 

DO 50 IGAUS-t rNCAUS 
DO SO JGAUS-UHfeAl^ 
LGAUS-LGAUS-i«l 
XEOI V-POSGt t I|AUS> 
YEOIV-PO^Gt < JfeAUBT 







c 

C mE «^U8S fOSITIONS TO CALCULATE LOCAL VALUES 

C 

CALL SKA**EAt 0ESIV,SHAPE,XEQIV, YEQIV) 

C 

C SET Uf* JACOBIAN HATRIX AND INVERSE 

C 

CALL D JACOB ( DER I V , DET J6 , D J ACK , D JAC I , NNODL , 

I IELE«*NELEW) 

C 

C CALCULATE CLOBAL DERIVATIVES FOR LINEAR FUNCTIONS 

C 

DO 30 tDinE»1,E 
DO 30 IN0DL«1,4 
CARTL( lOIHE, 1NODL>»0.0 
DO 10 JD1HE«I,E 

CARTLC IDIN£> INOOL >»CARTLt IDIHE, INDDLI’I’OJACl (IDINE, JDINE)* 

I DERlV<JDinE,INODL> 

10 CONTINUE 

30 CONTINUE 

DUOXt«0.0 
DU0Y1«O.O 
OVDXIxO.O 
DVDY1«*0.0 
DO 32f IN0DL«t,4 

DUDXI-DUDXI+CARTLC 1 , 1N00L)*UN0DE< lELEW, INODL) 

DUDY 1 -DUDY1 +CARTLC S, INODL ) ♦UNODE ( lELER , INODL ) 

DVDX I »DVDX I +CARTLC 1 , INODL ) ♦VNODE ( lELEM, INODL ) 

DVDY1 ■DVDYI+CARTLCE, INODL l•VNODE( lELEH, INODL) 

3Bf CONTINUE 

SUH1 aCAUS)-DUDX! 

SUflE«LCAUS)«DVDX! 

8UH3<L«AUS)-DUDY! 

SU»4(LCAU8>-DVDYI 
$UKltLSAUS)»0.0 
DO IS INODL»t«A 

SUHI C LCAUS ) •SUHI ( LSAUS ) •(•RHONOD ( I£LEM« INODL > pSHAPE ( INODL ) 

IS CONTINUE 

8UHy<LCAUS)»0.0 
DO 17 lN0DL«t,4 

SUNY ( LC AU8 ) »SUHY ( LCAU8 1’l'EOSS ( lELEH « INODL ) ♦SHAPE C INODL ) 

17 CONTINUE 

ARE AW < LC AU8 ) »DET JBA WE I CL C IC AUS 3 ♦WE 1 CL ( JCAUS ) 

C 

C CALCULATION OF PLASTIC WORK FUNCTION 

C 

CALL PL AST I C < LC AUS , SUH 1 , SUM2 , 8UH3 , SUH4 , 8UH Y , 8UH I , ARE AW , U ) 

50 CONTINUE 

SUH»<W(f )i‘Wt2)+W(3)+W(4)) 

RETURN 

END 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SU BROUT I NE PIRST I C i LCSUS 1 8 UH I , SUH2 , SUM3 , SUH4 , 8UHY , 8UH I , 

1 AREAW W) 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

IMPLICIT REAL • 8(A-H,0-Z) 

DIMENSION SUM1 t4),SUM2{4)»SUM3(4),SU«4t4), 
t syMI(4).A»EAWt4),W(4».8UBY<4) 

U>*<SUM1 {LCAUS)-8WM(LCAUS))**E 
M2* < SyM2 ( LCAyS )+8yM3 < LC AOS ) ) **2 
W3*8UH1 ( LCAOS }+SWt4(UCAtfS> 

W4»SyMt(LCAUS) 




rji r» 


IIS»$UK4a€AUS^ 

y^-syHKtCAus) 

m» «fi.o/9.o)**<yi ♦y4*»E+y5#^2**- 

t ,0-W6)**»O.5t4 
Wf • i W3/ i a , 4»*Me > )♦♦£ 

SOT I *43 8UHY i tCAUS )+0.0!)<«*0.32) 

y { ty Aus ) "* < y6«>*a « 5 > ( y7'«-u9 ) * 

I AaEAyCLCAU8)*«>$a7t 
RETURN 
ENU 

«iii##««#«i«#«»iit»i««*«t«#««»«t»*tit»*««*«**«#»««*«»*««*i*iii«ii 



AiJovo,' 
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